o 
o 

(N 
cn 



Science with a lunar low-frequency array: 
From the dark ages of the Universe to nearby exoplanets 

Sebastian Jester 

Max-Planck-Institut fiir Astronomie, Kdnigstuhl 17, 69117 Heidelberg, Germany 

Heino Falcke 

Department of Astrophysics, Institute for Mathematics, Astrophysics and Particle Physics, Radboud University, P.O. Box 9010, 6500 GL Nijmegen, The 

Netherlands 

ASTRON, Oiide Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands 



Abstract 

(-H Low-frequency radio astronomy is limited by severe ionospheric distortions below 50 MHz and complete reflection of radio waves below 10-30 
MHz. Shielding of man-made interference from long-range radio broadcasts, strong natural radio emission from the Earth's aurora, and the 
I need for setting up a large distributed antenna array make the lunar far side a supreme location for a low-frequency radio array. A number 
O of new scientific drivers for such an array, such as the study of the dark ages and epoch of reionization, exoplanets, and ultra-high energy 
cosmic rays, have emerged and need to be studied in greater detail. Here we review the scientific potential and requirements of these new 
^ scientific drivers and discuss the constraints for various lunar surface arrays. In particular we describe observability constraints imposed by 
i__ifhe interstellar and interplanetary medium, calculate the achievable resolution, sensitivity, and confusion limit of a dipole array using general 
scaling laws, and apply them to various scientific questions. 

Of particular interest for a lunar array are studies of the earliest phase of the universe which are not easily accessible by other means. 
^ These are the epoch of reionization at redshifts z=6-20, during which the very first stars and galaxies ionized most of the originally neutral 
intergalactic hydrogen, and the dark ages prior to that. 

' For example, a global 21-cm wave absorption signature from primordial hydrogen in the dark ages at z = 30-50 could in principle be detected 
by a single dipole in an eternally dark crater on the moon, but foreground subtraction would be extremely difficult. Obtaining a high-quality 
. power spectrum of density fluctuations in the epoch of reionization at z = 6-30, providing a wealth of cosmological data, would require about 
10^-10^ antenna elements on the moon, which appears not unreasonable in the long term. Moreover, baryonic acoustic oscillations in the dark 
ages at z=30-50 could similarly be detected, thereby providing pristine cosmological information, e.g., on the inflationary phase of the universe, 
^ij ' With a large array also exoplanet magnetospheres could be detected through Jupiter-like coherent bursts. Smaller arrays of order 10^ antennas 
over ~ 100 km, which could already be erected robotically by a single mission with current technology and launchers, could tackle surveys of 
^ steep-spectrum large-scale radio structures from galaxy clusters and radio galaxies. Also, at very low frequencies the structure of the interstellar 
ftiedium can be studied tomographically. Moreover, radio emission from neutrino interactions within the moon can potentially be used to create 
^ . ^ neutrino detector with a volume of several cubic kilometers. An ultra-high energy cosmic ray detector with thousands of square kilometre 

area for cosmic ray energies > 10^" eV could in principle be realized with some hundred antennas. 
" " In any case, pathfinder arrays are needed to test the feasibility of these experiments in the not too distant future. Lunar low-frequency arrays 
are thus a timely option to consider, offering the potential for significant new insights into a wide range of today's crucial scientific topics. 
This would open up one of the last unexplored frequency domains in the electromagnetic spectrum. 
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1. Introduction 
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Low-frequency radio astronomy is currently experiencing an 
impressive revival with a number major new facilities in opera- 
tion or under construction, such as the Giant Metrewave Radio 
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Telescope (G MRT.ISwarup et a Ull99l HIi Low-Frequ ency AR- 
ray (LOFAR, iFalcke et alf liooe; Ro ettgering et all 2006F , 



Long-Wavelength Array (LWA, Kassi m et al. 



Murchison Widefield Array (MWA, f ormerl y known as the 
Mileu ra Widefield Array, iMoralesL |2005[ [Morales et all 



2006l{'*L the 21 Centimeter Array (21CM4£J; formerly called 
Primeval Structure Telescope, PASlf^ Peterson et aL, 2005) . 
and P recision Array to Pro be the Epoch of Reionization (PA- 
PER, iBradley et al.. '2005"). Eventually the Square-Kilometre 
An-ay (SKA, [Schilizzi, 2004; Carilli and RawHngs, 200413 
may extend the expected collecting area even further. Main 
science drivers of these telescopes are the so-called epoch 
of io nization (e.g.. Shaver and de Bruyn, 2000; Carilli et al., 
2004 and references therein), large scale surveys , transient 
and variable source monitoring ('Fende r et al. , 2006h . observa- 
tions of our own sun (Kasper, 2006) and so lar system (|Zarka, 
2000h . a s well as exoplanets jZarkal 120071) and astroparticle 
physics dFalcke and Gorhaml l2003h . This rapid development 
together with the lively discussion of a possible return to the 
lunar surface by various space agencies, has raised again the 
long-standing question about the potential of low-frequency 
astronomy from space (e.g., Lazio et al., 2007). 

The currently planned ground-based telescopes will provide 
a serious advance in radio astronomy and extend the accessible 
frequencies to the widest range possible from the ground. The 
Earth's turbulent ionosphere gives rise to "radio seeing", mak- 
ing ground-based radio observations of the sky become diffi- 
cult at V < 100 MHz. At even longer wavelengths below about 
10-30 MHz one encounters the ionospheric cut-oflt where ra- 
dio waves are reflected, permitting long-distance short wave 
transmission around the Earth, but prohibiting observations of 
the sky. Observing just above the cut-off, i.e., between ~10- 
50 MHz requires especially favourable geomagnetic and iono- 
spheric conditions to obtain any decent images. The range be- 
low the cutoff is only readily observable from space. Hence, the 
dominant "low-frequency/long-wavelength" regime for which 
ground-based telescopes are being designed is at frequencies 
> 30 MHz and wavelengths < 10 m. We will therefore hereafter 
refer to wavelengths above 10 m as the ultra-long wavelength 
(ULW) range. This wavelength-based definition is preferable 
to a frequency-based one to avoid confusion with the technical 
very low frequency (VLF) designation that extends only from 
3-30 kHz; the range 300 kHz-30 MHz is in engineering lan- 
guage named the medium and high frequency regime (MF & 
HE), which is clearly not appropriate in a modern astronomical 
context. 

The best resolution achieved so far in the ULW range is 
on the scale of around 5 degrees, but more typically of or- 
der tens of degrees. This compares rather unfavourably to the 
milU-arcsecond resolution that can be routinely obtained in very 



' |http : //www ■ gmrt . ncra . tif r . res ■ in| 
^ http://www.lofar.org 
^ http://lwa.urmi.edu 

[http : //www . haystack . mi t ■ edu/ast/arrays/mwa/| 
^ 'http : //2 Icma . bao . ac . cnT* 

* http : / /web . phy s . emu . edu/~p ast/overview . html| 
^ |http : //www ■ skatelescope . org] 



long baseline interferometry (VLBI) at higher radio frequen- 
cies. Hence, the low-frequency Universe is the worst-charted 
part of the radio spectrum, and perhaps even of the entire elec- 
tromagnetic spectrum. By today, only two kinds of maps of 
the sky have been made at frequencies below 30 MHz. The 
first are maps of a part of the southe rn sky near the Galactic 
center such as tho s e obta ined by Cane and Whitham ( 1 9JJ), 
EUis and Mendillol (Il987l) and lCane and Ericksonl(l200lh from 



Tasmania. These have angular resolutions ranging from 5 to 30 
degrees. T he second kind are the map s obtained by the RAE- 
2 satellite (iNovaco and Browiil 19781) with angular resolution 
of 30 degrees or worse. None of these maps show individual 
sources other than diffusive synchrotron emission of the Galaxy, 
nor do they cover the entire sky. 

To improve this situation and to overcome these limita- 
tions, space-based low-frequency telescopes ar e requi r ed for 



all observations be l ow the ionospheric cutoff dWeilei , 1987 



Weiler et all Il988l: iKassim and Weileii Il990l) . This is also 



true for a significant part of the seeing-afFected frequency 
range above the cutoff frequency where high-resolution and 
high-dynamic range observations are required, such as imag- 
ing of 21 -cm emission of n eutral hydrogen in the very early 
Universe dCarilh et alll2007l) . 

So far, there have only been two space missions whose pri- 
mary purpose was low-frequency radio astronomy: the Radio 
Astronomy Explorers (RAE) 1 and 2. The design and results 
of these missions have been reviewed by lKaiseii (ll987h . lKaiser 
(1990) summarized the science goals of the RAE missions as 
four-fold, namely characterizing the spatial and spectral struc- 
ture of Galactic radio emission, observing bursts from the Sun 
and Jupiter, and attempting to detect discrete Galactic and ex- 
tragalactic sources. It was a surprising finding of RAE-1 that 
the Earth itself is a strong emitter of low-frequency bursts, the 
so-called Auroral Kilometric Radiation (AKR) generated by 
solar-wind interactions with the Earth's magnetic field. This 
emission is so strong that RAE-2 was placed in a lunar instead 
of the originally planned terrestrial orbit to provide shielding 
from this natural terrestrial interference (Fig. [T] and § l2.3l l. Of 
the original science goals, only the solar and Jupiter burst stud- 
ies could be completed; in addition, terrestrial bursts were stud- 
ied. The structure of the Galaxy's emission was only seen at 
very low spatial resolution (with beam sizes between 37° x 61° 
and 220° x 160°) and fairly low signal-to-noise ratio. Due to the 
very limited angular resolution and the large power of AKR, it 
proved impossible to image any discrete sources directly. Thus, 
two of the initial scientific questions of the RAE missions re- 
main unanswered to date. 

Low-frequency radio imaging telescopes at these frequencies 
are nowadays realized as digitally phased arrays of thousands 
of simple dipole antennas that are spread over a large area. This 
puts a serious constraint on any realization consisting of free- 
flying antennas in space and very likely favours lunar surface- 
based concepts in the long term. Moreover the moon provides 
unique shielding against disturbing radio interference from sun 
and Earth and within the exploration efforts of space agencies 
lunar missions appear to be a realistic possibility. Hence we 
will focus here on a lunar surface telescope concept consisting 
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Fig. 1. Shielding of terrestrial radio interference by the moon, as observed by the RAE-2 satellite. Even at a distance of 300 ,000 km from the Earth, the additional 
shielding by the Moon provides 1-3 orders of magnitude (10-30 dB) of additional interference suppression. Taken from [Alexander et alj ) 19751 © ESO). 



of dipoles as the most basic antenna unit. Nonetheless, most 
of our conclusions are general enough that they can be of use 
also for other kinds of implementations. 

Over the past four decades several conferences and 
study groups formed by space agencie s have produced re- 
ports on lunar radio astronomy (■Weiie?, '1987'; 'weiler et al.', 
.1988; Kassim and Weiler, 1990; Nieuwenhuizenet_aL, 1992; 
Landecke r et all Il99ll. Il99 2; Daint vt il992l; iFoingj |l994, 
1 QQA- Belv e* fiQQTl 'vxT^n 1 007.- iTrvnoc ot oi I 'iQQsl nnnn- 

et all 



1996; Bel 
Stone 



200C 



19971; 'Woan, 19 97t Jones et al.Vl998V2000; 
Takahashi, 2003r iOberoi and Pincon, 2005; 



Kassim et al 1 l2005al) . Starting from these earlier studies we 



will here summarize the prospects and fundamental limitations 
of ULW observations and discuss their scientific potential in 
light of our current knowledge. Recent scientific advances in 
many areas and in particular in cosmology lets low-frequency 
astronomy shine in a new light and makes this a timely topic. 

The plan of the paper is as follows. We describe observabil- 
ity constraints for ULW observations from the Earth and the 
Moon in ^ including a description of the Galactic foreground 
emission, and formulae for the estimation of array sizes, sen- 
sitivities etc. We give science questions to be addressed by a 
lunar low-frequency array in ^and summarize our findings in 

m 



2. Observability constraints for ULW observations 

2.1. Earth Ionosphere 

Radio waves are scattered by free electrons, such as those 
occurring in a space plasma or the Earth's ionosphere. The scat- 



tering strength increases towards lower frequencies. At frequen- 
cies below the plasma frequency, which increases with elec- 
tron density, radio waves cannot propagate at all. The plasma 
frequency of the Earth's ionosphere is typically near 10 MHz 
by day and near 5 MHz by night, implying that radio obser- 
vations near and below these frequencies are impossible from 
the Earth (the plasma frequencies can reach somewhat lower 
values near the Earth's magnetic poles). At frequencies up to a 
few 100 MHz, the ionosphere causes celestial radio sources to 
suffer significantly from angular displacements, angular broad- 
ening, and intensity fluctuations (scintillation). These features 
are akin to seeing, the twinkling of stars visible to the eye, 
which is caused by density fluctuations of turbulence layers 
in the Earth's atmosphere. Because of the scattering, the best 
available ground-based maps of the low-frequency sky have 
resolutio ns ranging from 30 degrees at 2.1 MHz to 5 degrees at 



resolutio ns ranging irom i\j aegrees at z.i MMz to 3 aegrees at 
10 MHz dElhs and Mendilloll 19871: ICane and Ericksonll200lh 
and very poor dynamic range. 

2.2. Lunar Surface and Ionosphere 



Even though the moon does not have a gaseous atmosphere, 
instruments carried by the Apollo 14 and Luna 19 and 22 
missions showed that the moon's surface has a photoelec- 
tron sheath, i.e., a weak ionosphere. On the day side, the 
inferred electron densities lay in the range 500-10,000/cm^, 
with plasma frequencies between 0.2 and 1 MHz (see 



Reasoner and Burke 1972, Benson et al. 1975 and references 



therein). Electron densities on the night side are expected to 
be much lower, since the moon's surface potential becomes 
negative. However, the electron density has not been mapped 
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yet in detail as function of lunar position and altitude above 
the surface. This is clearly a necessity, given that electron 
densities at the high end of the inferred ranges would mean 
that a lunar observatory only gains a factor of 10 in frequency 
during lunar day compared to what is observable below the 
Earth's ionosphere. The current generation of moon missions, 
in particular the Japanese KAGUYA/SELENE mission, should 
clarify this issue. 

Secondly, low-frequency radio waves can propagate into the 
lunar regolith to depths of 1-100 wavelengths. Subsurface dis- 
continuities in the electrical properties of the lunar regolith 
could therefore lead to reflections and hence stray signals be- 
ing scattered into the beam of lunar antennas. Thus, the sub- 
surface electrical properties of the prospective sites for a lunar 
radio telescope need to be explored before construction of a 
very large array (compare iTakahashL ,2003,) . 



2.3. Man-made and natural Interference 

The low-frequency radio spectrum is occupied to a large frac- 
tion by terrestrial broadcasts, in particular FM and longer-wave 
radio and television broadcasts and military and civil commu- 
nications. This radio frequency interference (RFI) is a severe 
problem for all kinds of ground-based radio telescopes. These 
signals block out cosmic signals directly at the correspond- 
ing wavelengths, but they also lead to an increased noise level 
for observations at other frequencies through intermodulation 
products in the receiving system, through leakage, or if the af- 
fected frequencies are not filtered out completely. 

At the longest wavelengths, these signals are reflected by 
the ionosphere. This is of course the very reason why these 
wavelengths are used for world-wide radio broadcasts, but it 
also means that a terrestrial ULW radio telescope is sensitive 
not just to local interference, but also to interference from all 
parts of the world, irrespective of location. 

When the first dedicated radio astronomy satellite, RAE-1, 
was designed, it was expected that placing the satellite in a 
low orbit above the ionosphere would lead to a sufficient sup- 
pression of interference. However, it turned out that there were 
very strong radio bursts from the Earth itself, generated by the 
interaction between solar-wind particles and the Earth's mag- 
netic field, in particular in the frequency range 150-300 kHz. 
At higher frequencies, man-made and lightning interference are 
also still detectable above the ionosphere. Therefore, the sec- 
ond mission, RAE-2, was redesigned and this satellite was put 
into an orbit around the Moon. Figure [T| shows how well ter- 
restrial radio interference is shielded by the Moon. 

Naturally, also the sun and Jupiter are strong sources of ULW 
radio emission - especially during outbursts - and again here 
the moon can act as a perfect shield for a large fraction of the 
time. 



attenuation around Malapert Mountain 



0.5 MHz 
1.0 MHz 



Takahashil (l2003h has modelled the reflection and absorp- 




tion of low-frequency radio waves by the Moon, including 
the shielding by mountains. He finds that a large fraction of 
the incident energy is reflected by the lunar surface (55% at 
0.5 MHz), and an eff'ective skin depth of 2.2km(v/MHz) " 



-210 -180 -150 -120 -90 -60 -30 
distance from South Pole (km) 

Fig. 2. Shielding of terrestrial radio interference by Malapert mountain, whose 
peak is at x 120 km from the south pole. This shielding is in addition to the 
20-30 dB at tenuation applyin g to all interference in the south polar region. 
Taken from IXakahashl j2003l Fig. 4.9) with kind permission of the author. 

Thus, at 0.5 MHz, shielding is already substantial at a depth 
of a few kilometers and becomes better at higher frequencies. 
Figure |2] shows the attenuation along a line joining the lunar 
south pole (a possible location for a future lunar base) and the 
nearby Malapert mountain; Malapert may provide up to 70 dB 
(a factor of 10^) of shielding at 0.5 MHz, while the simulations 
indicate a shielding by a factor of at least 10** at all frequencies 
above 10 kHz for the lunar far side, i.e., substantially stronger 
shielding. 



2.4. The interplanetary and interstellar media 

A fundamental constraint for all very low-frequency radio 
astronomy is the presence of free electrons in the interplanetary 
medium (IPM) within the solar system, the solar wind, and in 
the interstellar medium (ISM) filling our Galaxy. The plasma 
frequency of the IPM depends on the distance from the sun 
and the solar cycle; at the Earth, it lies in the range 20-30 kHz, 
and the scaling is roughly linear with sun-Earth distance. Ob- 
servations outside the solar system are restricted to higher fre- 
quencies. The ISM plasma frequency is of order 2 kHz for the 
average electron density in the Galactic plane of 0.025 cm~^ 
dPeterson and Webben, l2002h and therefore poses no additional 
constraints beyond those imposed by the IPM. The ISM plasma 
frequency is significantly lower than the ionospheric plasma 
frequency and hence there is still a wide range of accessible 
parameter space to be explored. 

However, in addition to inhibiting the propagation of elec- 
tromagnetic waves below the plasma frequency, there are 
several more consequences of the presence of free electrons: 
foreground emission, free-free absorption, and scattering of 
waves in turbulent regions. As usual all these effects can be 
a curse, when studying extragalactic sources, and a blessing, 
when studying the local ISM. Moreover, there are fundamental 
constraints to observations of the radio sky with interferom- 
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Fig. 3. Spectrum of Galactic synchrotron emission, sliowing specific intensity 
(solid line) and surface brightness temperature (dashed line). The surface 
brightness temperature do minates the system t e mpera ture at the frequencies 
shown. Figure taken from lOberoi and PinconI i2003h with kind permission 
of the authors. 

eters, which arise from the need to match the resolution to 
the density of sources at the instrument's sensitivity to avoid 
confusion. We consider each of these effects in turn. 



2.4.1. Galactic synchrotron emission 

Synchrotron emission from electrons moving in the Galac- 
tic magnetic field constitutes the dominant foreground (an 
overview about observations of the Galactic emission that 
have been made to date is given by DwarakanathI 2000l and 
Keshet et al.l 120041) . The distributions of electrons and mag- 
netic fields are not yet known in full detail and subject of 
active study. The component of the ISM that contains the free 
electrons responsible for the scattering is called the Warm 
Ionized Medium (WIM). It broadly follows the distribution of 
gas in the Galaxy, i.e., it is concentrated in the plane of the 
Galaxy. The brightness temperature of the emission rises from 
about 10"* K at 30 MHz to about 2.6 x 10^ K at around 1 MHz, 
and then increases less rapidly (and the specific intensity de- 
creases) because free electrons absorb low-frequency radiation 
via free-free absorption (see Figure [3] and following section). 
For reference, the conversion between flux density (Jansky, 
JyI3 and brightness temperature (in Kelvin, K) is given by 
Sv - ByQ. - Ikc^^v^T^Q. and hence we have for an extended 



5v = 0.93kJy 



Q 



lO^K deg^ 



1 10 MHz/ 



(1) 



where By is the surface brightness, is the surface brightness 
temperature, 5 V is the flux density, Q. is the source's solid angle, 
and k is the Boltzmann constant 1.38 x 10"^^ J/K. 

Because of the absorption, the surface brightness distribu- 
tion of the Galaxy varies with frequency: at high frequencies, 
the Galactic plane is seen in emission, with HII regions being 



superimposed in absorption, while the poles are more or less 
transparent; at frequencies around 2-3 MHz, the surface bright- 
ness is constant in all directions, so that the sky appears equally 
foggy everywhere; at lower frequencies, the plane is seen in 
absorption, so that the poles will appear warmer than the plane. 

The main effect of the Galactic synchrotron emission at low 
frequencies is that the sky temperature is so high that it domi- 
nates the system temperature and therefore limits the sensitiv- 
ity. Figure |4] shows the expected sensitivity for a lunar array as 
function of the number of crossed dipoles, each with effective 
area 



Aeff = ^74, 



(2) 



where effective area is related to primary beam size Q, which 
is also the size of the antenna's field of view, by 



(3) 



We use the RMS sens itivity (Icr) of an array of antennas 
( [Thompson et al.ll200ll the origin of this expression is the sys- 
tem equivalent flux density, SEFD, which is the surface bright- 
ness corresponding to the system temperature) as given by 
Cohea (.2004) : 



o"RMS =58.9 mJy/beam X 



AeffVA^(A^-l)MFfexp AV 



(4) 



where Aetf is the antenna's effective area in m^, Av and fgxp are 
the bandwidth in MHz and observing time in hours, respec- 
tively, over which the signal is integrated, and is the number 
of array elements (for a single dipole, the term N(N - 1) needs 
to be replaced by 1). T^ys is the system temperature, i.e., the 
sky temperature (given in Fig. O and approximated here as 

163x10^ k(—^) 
\2MHz/ 

16.3X io^/:(— ^1 

\2MHz/ 



! sky 



-2.53 



at V > 2 MHz 



at V < 2 MHz. 



(5) 



For a dish antenna, A^ip is the number of intermediate frequency 
(IF) bands, i.e., 1 or 2 depending on whether only one or 
both polarizations are recorded simultaneously; for our crossed 
dipoles, the effective area already includes the factor of 2 for 
recording both polarizations simultaneously, so that A^ip = 1 is 
appropriate here. 

For later reference, we can also cast eqn.|4]in terms of surface 
brightness sensitivity: 



RMS - 



^max -' sys 



Aeff y/N{N - 1) MpfexpAv 



T 

sys 



(6) 



(7) 



where Z)n,ax is the maximum baseline of the array setting the 
maximum angular resolution and 



/ 



Aeff VA^(A^- 1) 



(8) 



IJy = 10 WHz- 



10- 



ergs 



Hz- 



is the filling factor of the array, i.e., roughly the fraction of the 
arrays's physical area that is filled with effective antenna area. 
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Fig. 4. RMS sensitivity (1-sigma noise level) of an array of crossed dipoles (antenna effective area = /!^/4), for 5% fractional bandwidth and 24h integration 
time, as function of the number of crossed dipoles in the array as given by eqns. @ and j6); see eq. ^9) for the maximum useful number of dipoles due 
to shadowing. The system temperature is assumed to be equal to the Galactic sky's brightness temperature in the direction of the Galactic poles (Fig.[3)- 
The left-hand axis gives the point-source sensitivity. The right-hand axis yields the surface brightness sensitivity in Kelvin for one resolution element after 
multiplying the given value by the square of the an'ay's maximum baseline D, where D is measured in meters. The RMS scales as the inverse of the square 
root of both bandwidth and integration time. 



There is a fundamental maximum to the useful number of 
dipoles within a given area: an array with resolution arcmin- 
utes can have at most 



= 4.7 X 10^ ( p 



(9) 



dipoles before the filling factor exceeds unity, i.e., the dipoles 
start becoming mutually coupled; the maximum baseline nec- 
essary to reach this resolution is given by 



3438/1 



(10) 



We consider a related issue in ^2. 5 1 below: the maximum num- 
ber of useful antennas before confusion noise (rather than ran- 
dom noise) limits the array performance. 



2.4.2. Free-free absorption 

The magnitude of the free-free (or thermal bremsstrahlung) 
absorption coefficient Ky depends on both electron density and 
temperature as 



Ky = 1.78 ^ff 



■ill 



(11) 



where is the density in cm and T,, the temperature in K 
of the electron distribution, v is the observing frequency in Hz, 
and gs is the free-free Gaunt factor 



gs = 10.6+ 1.6 log r,- 1.26 log V 



(12) 



(IPeterson and Webben l2002r) . For the parameters of the 



interstellar medium n. 



0.025 cm"^' and 



7000 K 

(iPeterson and Webbeit |2002|) . the distance at which the ISM 
optical density satisfies t = 1 is about 1 .5 kiloparsec (kpc) at 
1 MHz; this distance scales as the square of the frequency. For 
comparison, the thickness of the Galactic disc, in which the 
free electrons are found, is about 1 kpc, and the diameter of 
this disc is about 25 kpc. Thus, extragalactic observations are 
impossible below about 3 MHz unless there are low-density 
patches. Indeed, the Galactic synchrotron spectrum (Fig. |3]l 
turns over at frequencies near 3 MHz due to absorption (this 
turnover is one of the ways in which the Galactic electron 
density and temperature are constrained observationally). Vari- 
ations of the electron temperature and density in the ISM will 
already cause difficulties in interpreting data in the range 3- 
7 MHz, since variations in the surface brightness arising from 
fluctuations in the foreground emission and absorption have to 
be disentangled from structures in the sources. For some kinds 
of sources, this may be done by comparing the source spectra 
as well as their morphology, but for extragalactic synchrotron 
sources with similar spectra to the Galactic synchrotron emis- 
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sion, this task will be considerably more difficult the more 
extended they are. 

On the other hand, the absorption enables a tomography of 
the ISM, since observing at different frequencies picks out ab- 
sorbing structures at different distances. This is described in 
section 13.3.11 

2.4.3. Brightness temperature and dynamic range 

Below about 100 MHz, the brightness temperature of the 
Galactic emission is much higher than that of any noise source 
in the telescope itself and therefore sets the system temperature, 
Tsys. Because of the structure of Galactic emission and absorp- 
tion, Tsys at any given frequency will vary strongly with direc- 
tion (for example, below 30 MHz, HII regions w i ll be seen in 
absorption against the Galactic plane). Ericksonl ( 119991) noted 
that this can lead to dynamic range and calibration problems 
already when observing with the VLA's steerable antennas at 
74 MHz; these difficulties are likely to be exacerbated by the 
use of lower-frequency dipoles with larger primary beams. The 
operation of LOFAR will provide experience and expertise in 
the calibration of an interferometer in the presence of strongly 
variable backgrounds. For some applications a dynamic range 
of > 10^ is sought. Since the field of view of dipole antennas is 
practically the entire sky, the necessary dynamic range can be 
estimated by comparing the flux densities of the brightest ra- 
dio sources in the sky, the supernova remnant Cas A (diameter 

6', flux density 65 kJy) and the radio galaxy Cyg A (diam- 
eter ss 2', flux density 31.7kJy; both flux densities are taken 
from Tab. 2 of iBaars et al.l Il977h . to the desired flux limit. For 
extragalactic sources, the latter is the confusion limit (see eq. 
l22l below). At 15 MHz and with a resolution of 2', for exam- 
ple, the confusion limit is 0.14Jy, requiring a dynamic range 
of X! 2 X lO*" to image sources at the confusion limit against the 
glare of Cyg A. 

2.4.4. Scattering in the Interplanetary and Interstellar Media 
(IPM/ISM) 

Like all solar-system planets, the Earth is embedded in the 
interplanetary medium (IPM): the solar wind, an extension of 
the solar corona. The properties of the solar wind vary with 
the 11 -year solar activity cycle. The solar system itself is em- 
bedded in the ISM, as already discussed above. Electron den- 
sity fluctuations both in the ISM and the IPM cause scattering 
of radio waves, leading to the phenomena of diffractive and 
refractive scintillation, scatter broadening, temporal broaden- 
ing of pulses, and Faraday depolarization of linearly polarized 
sources. The angular scattering and scintillation are akin to 
"seeing". This is the blurring of optical images in ground-based 
telescopes caused by density fluctuations in turbulent layers 
of the Earth's atmosphere and the main motivation for optical 
space telescopes and adaptive optics. 



2.4.5. Angular scattering 

At frequencies below 30 MH z, the scattering is in th e 
"strong" regime for both media dCohen and Cronvnl 1 19741) . 
which means that there are no intensity scintillations, and the 



scattering broadens the intrinsic source size into a scattering 
angle that scales roughly as A-. There is solar-wind scattering 
at all elongations from the sun, but the scattering is smallest in 
the anti-solar direction and increases for lines of sight passing 
closer to the sun. Similarly, the ISM is completely surround- 
ing the Earth so that there is some scattering in all directions, 
but because of structures in ffie ISM, the scattering is stronger 
for lines of sight passing through the plane of the Galaxy and 
closer to the Galactic centre. 

The expected order of magnitude of ffie scattering for both 
media can be extrapolated from measurements at higher fre- 
quencies. Empirically, the ISM broadening of the image of an 
extragalactic source is roughly given by 



&1SM 



30' 



(v/MHz)2-2 Vsin^ 



(13) 



where b is the Galactic latitude (iCohen and Cronyni Il974h . 
Estimates of the IPM scattering at 1 MHz range from a fac tor 
3 smaller than the ISM scattering dRickett and Colesl I2OO0I) to 
a factor 5 greater (Woan, 200Ci). Figure |5] shows the achievable 
resolution with the limits imposed by IPM and ISM scattering, 
as function of array size, using the worst-case IPM scattering 
given by Woan (2000,): 



100' 



(v/MHz) 



2' 



(14) 



and a slightly different version of eq. [13] for ISM scatter- 
ing, 19-i.sM - 22'(^ v/MHz)^-^. The best-case IPM scattering 
dRickett and Colesl |2000.) would be a factor of 15 smaller than 
given by eq. dT4l i. 

The scaling i? oc arises from a thin-screen approximation 
for the scattering medium. For a screen of thickness L with 
overdensities of characteristic radius a and excess density AA^^, 
the scattering angle is given by i^thin = r^A^ANe'^/iL/a)/!?!, 
where = 2.82 x 10"'^ m is the classical electron radius. A 
more realistic model of the ISM/IPM accounts for their turbu- 
lent nature, yielding the v^^ scaling used above. The strength 
of the scattering is set by the scattering measure, a weighted 
line-of-sight integral of the turb ulen ce power-spectrum coef- 
ficients (see Cordes et al. , 1991 and Cordes and Laziol 2002, 
who also give expressions relating ffie scattering measure to 
observables, as well as the emission measure and the pulsar 
dispersion measure; 'Cordes and Laziol 2002 furthermore re- 
port a model for the distribution of free electrons in the Galaxy 
and give a computer programme ffiat allows computation of the 
scattering measure, angular broadening, etc.). 

2.4.6. Temporal broadening 

Together with angular broadening, electron density fluctu- 
ations due to ISM and IPM turbulence also cause temporal 
broadening of radio pulses. The temporal broadening arises be- 
cause the angular broadening allows the emission to reach the 
observer via multiple optical paths with different lengths. The 
magnitude of temporal broadening Arb is just the geometric 
delay between the scattered and unscattered signal paths and 
therefore related to the magnitude of angular broadening by 
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Fig. 5. Resolution of a radio array as function of baseline. The ISM and IPM scattering limits are also shown, assuming the scattering laws given by IWoarJ 
jiooo). The ISM scattering size is in principle a firm limit to the resolution achievable for extragalactic observations, but could be smaller in some small part 
of the sky, in particular away from both the Galactic plane and the Galactic centre. The IPM scattering has not been measured well at low frequencies and 
could be a factor 15 smaller than given. The IPM scattering is given here for the anti-solar direction; it will increase when observing closer to the sun, and 
vary with the solar cycle. 



ATb = 



2c 



(15) 



in a small-angle approximation. Here z* is the reduced distance 
to the scattering screen and given by z* - (zz')/(z + z'), with z 
the observer-screen distance and z' the screen-source distance. 
For IPM scattering of Galactic sources and fo r all ex tragalactic 
observations, z' » z so that z* ~ Z- 'Woan| ( 2000h gives the 
following rough scalings for temporal broadening: 



ISM : Af = 6 yr (v/MHz) 
IPM:Af=0.1s(v/MHzr 



-4.4 



(16) 
(17) 



With these scaling laws, the IPM temporal broadening is ac- 
ceptable at least down to 1 MHz, but the ISM broadening is 
still 2 hours at 10 MHz, implying that no high-time resolu- 
tion work can be done for sources at extragalactic distances. 
For sources within the Galaxy, Af increases roughly with the 
square of the distance to t he source for fixed turbulence strength 
dCordes and Lazioll2002l eq. 9; in addition to the explicit linear 
dependence, there is a distance dependence in the line-of-sight 
integral for the scattering measure); this will be relevant for the 
observability of radio bursts from extrasolar planets 03.4.2l be- 
low). However, there are now indications (P. Zarka,pnv. comm.) 
that the pulse broadening may increase less rapidly with wave- 
length than expected from eq. ( fTSl i: this might allow the detec- 
tion of ULW pulses also from extragalactic sources after all. 



2.4.7. Faraday rotation and depolarization 

The plane of polarization of a linearly polarized beam under- 
goes the so-called Faraday rotation when it propagates through 
a region containing both free electrons and a magnetic field 
component parallel to the direction of propagation. Faraday ro- 
tation can be corrected using multi-frequency observations and 
the known frequency dependence of the effect. However, Fara- 
day rotation for low-frequency observations is dominated by 
free electrons from the solar wind in the immediate vicinity of 
Earth. This ma kes linear p olarization angles difficult to calibrate 
for all sources dWoan , 2000). since the RMS rotation angles can 
be very large and time- variable. On the positive side this may 
provide some measure of solar wind activity and distribution. 

A much more severe effect is the cellular depolarization aris- 
ing from multiple uncorrected Faraday rotation events in the 
ISM. The magnitude of the effect can be found by consider- 
ing the random walk of the polarization angle in a region with 
fluctuating electron density and magnetic fields, akin to the ran- 
dom walk of phase in a fluctuating electron density; assuming 
a constant ma gnetic field, the RMS phase variation is given by 
(IWoanLl2000h : 



Al/^; 



RMS 



2.6 X 10"'^ AN, 



A^a Bu 



LI a rad. 



(18) 



where again L is the path length through the scattering region, 
a is the size of an electron density fluctuation, AN, is the as- 
sociated electron overdensity, and By the magnetic flux density 
component along the line of sight (all in SI units). By inserting 



8 



into equation ( fTSl l a size a - 10"^ parsec for individual turbulent 
cells of the ISM, a typical electron density of 0.3 x 10^ m"^, and 
a typical magnetic field of 0.5 nT, the magnitude of the Faraday 
rotation at 1 MHz is about 1 radian in one such cell. The to- 
tal RMS Faraday rotation of low-frequency radiation travelling 
through the ISM is therefore very large, of order 1000 radians 
per parsec, and therefore leads on average to depolarization. 
Thus, no linear polarization is observable from sources beyond 
our own solar system at these very low frequencies. 

Circular polarization is not affected, but circularly polarized 
sources are extremely rare at higher frequencies, though cy- 
clotron maser bursts have been observed from certain types of 
stars (e.g. Hall inan et al.. 2007,) and is being found in an in- 
creasing number of active galactic nuclei dOabuzda et aUboOSl 
e.g.). At lower frequencies, cyclotron maser bursts are expected 
to be visible from extrasolar planets (see ^3.4.2l i. Also con- 
version fro m hnear to circular emission in a dense magnetized 
plasma (see iBeckert and Falcke , |2002[ and references therein 
for an explanation) is another process that can generate some 
level of circularly polarized emission across a wide frequency 
range. Hence, circular polarization capabilities - which can be 
measured by cross-correlating linearly polarized feeds - should 
be a more promising polarization signal to look for at low fre- 
quencies. 



2.5. Confusion and duration of a confusion-limited all-sky 
survey 

The first low-frequency radio surveys suffered badly from 
confusion: the presence of unresolved sources with individual 
flux densities below the detection limit leads to a constant floor 
in the noise level that is reached after a certain observing time. 
This can also lead to source posit ions being substantially in 
error (see, e.g., Kellermamj, [2003b . For low-frequency radio 
astronomy, a greater level of confusion may arise from fluctua- 
tions in the Galactic foreground emission rather than the pres- 
ence of unresolve d discrete extrag a lactic background sources. 

Among others, iDi Matteo et al. I (12004') have considered the 
impact of these fluctuations on the observability of the red- 
shifted 21 -cm power spectrum at frequencies near 100 MHz. A 
detailed simulation of the Galactic foreground emission based 
on an extrapolation of results from WMAP and LOFAR will 
provide useful guidance on the expected level of foreground 
fluctuations which need to be removed in order to detect ex- 
tragalactic sources. For the case of observing redshifted line 
emission, spectral information will help in isolating the signal 
from the foreground, and transient detection will be affected 
predominantly by the bright foreground emission; however, the 
observability of extragalactic continuum sources in the pres- 
ence of foreground fluctuations needs to be assessed in detail 
with suitable simulations. 

To estimate the confusion limit for extragalactic observa- 
tions, which can be don e down to 3 MHz, we use the approach 
and formulas of ICohenI ( 120041) . who estimates the anticipated 
confusion limit due to background sources by extrapolating 
from extragalactic source counts and the observed confusion 



limit at 74 MHz. As a conservative estimate based on compar- 
ing the RMS n oise level of 74 MHz observations in the VLA's 
C and D arravs. lCohenI (120041) requires a source density of less 
than one source per m - 12.9 beam areas for a survey to remain 
confusion-free. The beam area Qa of a Gaussian elliptical an- 
tenna beam with major and minor axes ctj and o",,, respectively, 
is given by 



Qa = 2n cr^cTy 
^ 1.13 )?l 



(19) 



where the approximation holds for a circular Gaussian beam 
(cTx - cTy - cr) of full-width at half-maximum (FWHM) & - 
V8 In 2 cr ~ 2.35 cr. The confusion limit S'conf is then defined 
by 

\.U'»^mN^{S,,,,f) = 1, 

where A^> (S ) is the number of sourc es with flux densities greater 
than S (equation 1 from Cohenj 2004). In the VLA Low- 
Frequency Sky Survey (VLSS) that is conducted at 74 MHz, 
the source counts at the faint end are given by 



N^iS) = 1.14deg 



-2 



(20) 



where is the number of sources per square degree, (3 - 1.3, 
and the flux density limit S is measured in Jy. To estimate N>(S) 
at still lower frequencies, we make the simplifying assumption 
that the source flux can be extrapolated in frequency with the 
spectrum v"°^. This is typical for optically thin synchrotron 
sources that are expected to dominate the extragalactic low- 
frequency population. Thence it follows that 



N>(S)^ ISOOdeg 



-2 



lOmJy 



"(— ) 

1 10 MHz/ 



-0.7 



(21) 



We ignore here the fact that many radio sources start showing 
a low-frequency turnover due to synchrotron self-absorption, 
and hence we likely overestimate the number of source and 
thus in turn also the confusion limit. With this assumption, the 
confusion limit in Jansky is given as function of resolution § by 

5eo„f('?,v) = (l2.9x 1.14x1.13 i^/rff"' (^^P 

= 16mJyx(./l')'-(^p. (22) 

The left panel of Figure|6]shows the confusion limit from equa- 
tion (I22I 1 as function of maximum baseline and frequency (see 
Fig.|5]for the resolution as function of those quantities). 

Since the confusion limit is a lower limit to the achievable 
noise, it implies an upper limit to the useful collecting area of 
an array — adding more antennas only decreases the time in 
which the confusion limit is reached, but not the overall array 
sensitivity. To compute the maximum useful number of anten- 
nas, we begin with the expression for the RMS sensitivity of an 
array of dipoles, eq. |4] As before, we assume that each array 
element is a pair of crossed dipoles with effective area given by 
eq. (|2]). Setting the confusion limit from eqn. l22l equal to the 



9 



Extragalactic, extrapolated from 74 MHz 

Tm| — I I I lllllj 1 I I lllllj 1 I llllllj 1 I llllllj 1 I I lllllj 1 I I lllllj — TTT 

3 



10' 
108 
10^ 
10* 
^ 103 

^ 100 

e 10 





1 


c 


.2 


0.1 


'm 




^0.01 


c 
o 


10-3 


o 






10-* 




10-5 




10-8 




10-' 




I- 1//MH 



iml I I I I I I 



10' 

108 

10^ 
10" 

103 

100 
10 

1 

0.1 
0.01 

10-3 

10-* 
10-5 

10-8 

10-' 



Extragalactic, extrapolated from 74 MHz 
10' 



6« 
m 
II 



108 ^ 



105 ^ 



10* r 



2 103 ^ 



(73 
ti-i 

c 
o 
o 



100 ^ 



10 r 



1 r 



TTn| — I 1 1 lllllj — I 1 1 lllllj — I 1 1 lllllj — I 1 1 lllllj — I I llllllj — I I llllllj — I 1 1 II 

, 300 I 




iz/MHz 

iml I I I I I I i_LLii 



10-3 0.01 0.1 1 10 100 1000 
Baseline [km] 



0.01 0.1 1 10 100 1000 10* 
Baseline [km] 



Fig. 6. Left: Confusion limit from extragalactic background sources as function of baseline. Right: Number of antennas reaching confusion limit at l-cr level 
in a 24 h integration with 5% bandwidth; the integration time to reach the same limit at 5cr significance would be 25 times longer, or about 1 month. 



1-cr RMS from eqn.|4]then yields the maximum useful number 
of antennas for a given integration time, which is plotted in the 
ri ght-hand pane l of Fig. |6]for an integration time of 24h. 
Perlev (1200 ih gives an expression for the "survey equation". 



the time needed to survey a solid angle of 1 sr to a given sen- 
sitivity. This is simply the time necessary for an array to reach 
the desired sensitivity (via eq. |4]i, divided by the field of view 
of each array element, i.e., the primary beam area. For the case 
of an array of crossed dipoles, with effective area of each 
element given by eq. |2] the time to survey one steradian to a 
1-cr point-source flux density s is given by 



ts 



urvey 



s 

44 d 



2 AvA^2^4 

-2 



Av 



T 

^ sys 

WK] \ 100 kHz 



\MHz/ 



(23) 



Inserting the high-frequency dependence of the sky temperature 
from eq. (|5]l, which dominates the system temperature, choos- 
ing the confusion limit from eq. (l22T i as limiting flux density, 
and including a telescope efficiency rji, this expression can be 
approximated as 



urvey 



3.3 d 



(-1 

1 100 J 



b 

oT 



(iMHz) 



-0.66 



-3.08 



(24) 



where b - Av/v is the fractional bandwidth, is the number of 
antennas, and fsurvey is the time necessary to survey one stera- 
dian of sky. This form of the survey equation is valid at fre- 
quencies V > 2 MHz and as long as the sky noise temperature 
over the system temperature; at frequencies below 2 MHz, the 
Galactic plane is nearly completely opaque, so that extragalac- 
tic sources cannot be observed in any case. An array of 300 



elements and maximum baseline 100 km, operating at 10 MHz 
(yielding a resolution of ^ 1') with b would thus be able 

to survey one steradian of sky (3282.8 deg^) to the confusion 
limit in 4/77^ days. Note that these equations use 1-cr sensitiv- 
ities; for higher significance levels, exposure times increase as 
the square of the desired significance level, e.g., for 5-cr sensi- 
tivities, the required time for an all-sky survey would be 3. 5/77^ 
years. Figure |5] shows the required exposure time and baseline 
as function of resolution for a confusion-limited all-sky survey 
with 5-cr significance at 3 MHz and 30 MHz. This is roughly 
the range that is accessible exclusively from above the Earth's 
ionosphere and not yet optically thick to free-free scattering. 

Equation ( |24] | gives an order-of-magnitude estimate. In prac- 
tice, details like (m, v) coverage will change the actual values. 
E.g., an array at a lunar pole will have more and more restricted 
(m, v) coverage for sources closer and closer to the lunar ce- 
lestial equator, losing the ability to resolve sources from each 
other in the direction perpendicular to the equator. 

As just mentioned, the presence of Galactic foreground emis- 
sion adds a new confusion component at frequencies below 
about 30 MHz. The accuracy with which extragalactic source 
observations can be done depends on the smallest-scale struc- 
tures that are present in the foreground emission. Observations 
of distant radio galaxies will be hampered by the fact that both, 
they and the Galactic foreground, produce synchrotron emis- 
sion, so that the spectral signature cannot be used for distin- 
guishing foreground and extragalactic emission, as is possible 
for 21 cm spectral line emission. As noted by Graham Woan, 
''There are clearly many contributions to the apparent flux den- 
sity of any single source, and much of the interpretatio n from , 
discrete sources will therefore be statistical in nature" ( Woanl 
|200Q|). However, at least at higher frequencies most of power 
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in the foreground is in the diffuse component and hence higher 
resolution clearly makes imaging more feasibly for any kind of 
survey work. 

2.6. Stability and Calib ratability 

The biggest challenge for a low-frequency phased array is 
the calibration. The main issue for terrestrial arrays is calibrat- 
ing ionospheric phase fluctuations. Due to the low-density at- 
mosphere and ionosphere of the moon - especially during the 
night - we do not consider this a major problem for a lunar ar- 
ray. This leaves calibration of the instrumental beam shape and 
band pass (frequency dependent gain of the antennas and elec- 
tronics) as the main concerns. For a phased array consisting of 
simple dipoles where all beam forming and imaging happens 
digitally, the overall beam pattern can in principle be calculated 
in a straightforward manner However, for terrestrial applica- 
tions relative gain changes in the analog parts of the dipole an- 
tennas need to be traced and calibrated out. For the low-band 
antennas of LOFAR the main antenna gain variations are ex- 
pected to come from temperature fluctuations (affecting elec- 
tronics and damping in the wires), rain (changing the electrical 
properties of ground and wires), and wind (leading to vibra- 
tions of the antennas). Actual variations of prototype hardware 
in the field at 40-80 MHz have been found to be of order 20% 



(iNehls et al.l l2007h . Since neither rain nor wind are to be ex- 



pected on the moon, temperature variations will be the largest 
source of gain variations. Lunar surface temperatures can vary 
from -233° C to 132° C, well in excess of terrestrial variations. 



Fortunately this change in temperature happens over the course 
of a lunar day, i.e., 4 weeks and should be relatively smooth 
and predictable over the array. Also, there is no variable cloud 
cover that can lead to random solar irradiation. The most sta- 
ble conditions on the moon would be achieved in an eternally 
dark spot in a polar crater, where one expects very little gain 
variations at all. 

The final source of uncertainty in the overall calibration 
would be the relative location of the antennas (if not known 
from the deployment procedure). However, interferometric self- 
calibration procedures can solve for antenna positions with high 
accuracy, especially if the sky is dominated by a single (known) 
point source. This is the case during Jupiter and solar bursts 
and has been used effectively alre ady to calibrate LO FAR pro- 



totype stations (such as LOPES; Falcke et al.L l2005h . 



Hence, in contrast to arrays on Earth, calibratability does not 
seem to be the main driver for design and layout of the array, as 
long as a proper temperature control of the dipoles is achieved. 

3. Science requirements 

This section describes science questions to be addressed by a 
lunar array. As far as possible, the science requirements (resolu- 
tion, sensitivity) have been translated into requirements on the 
array parameters (maximum baselines, number of antennas). 
However, there are additional requirements with regards to the 
two-dimensional distribution of antennas, the properties of the 
correlator, and the methods used to calibrate the interferomet- 
ric data and process them into images of the sky. These have 
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not been considered in detail and are likely to modify some of 
the array parameters presented here. Some of thes e additional 
aspects have been discussed by ' WoanI (Il997h and iJones et al.l 
(ll998), and other previous studies for lunar or space-based low- 
frequency arrays (see citations at end of section [U. 

3.1. Cosmology with HI line emission 

A fundamental question of current cosmological research is 
on the n ature of structure fo rmation in the Universe (for a re- 
view, see Ciardi and Ferraral 2005) : how is the observed distri- 
bution of visible matter created from the initial conditions just 
after the big bang, when matter and radiation were distributed 
extremely smoothly, with density variations of just one part in 
100,000? The CMB radiation was emitted at z a; 1200, about 
400,000 years after the Big Bang, when the Universe had cooled 
off sufficiently for electrons and protons to recombine into neu- 
tral hydrogen atoms (Epoch of Recombination), allowing the 
background radiation to move freely without being scattered by 
the electrons and protons. At the same time, however, the Uni- 
verse became opaque to visible light, because neutral hydrogen 
atoms absorb visible and infrared photons and re-emit them in 
random directions. Moreover, there were no sources of light in 
the Universe yet: the hydrogen and helium that were created in 
the big bang first had to cool in order to be able to clump to- 
gether and form stars and galaxies. Hence, this era is call ed the 
"cosmic dark ages" in the redshift range z = 30 - 1000 (iReesL 
[19991). 

Things only changed after the first stars, galaxies and active 
black holes had formed and emitted enough UV and X-ray 
photons to reionize the neutral hydrogen, allowing all radiation 
to pass freely. The time when this happened is called the Epoch 
of Reionization (EoR) and is believed to have occurred around 
z ~ 11, about 400 million years after the Big Bang, though it is 
at present not known whether the reionization happened more or 
less instantaneously, similar to a global phase transition, or was 
more or less spread out in time, depending on local conditions. 

Throughout all these epochs hydrogen played a major role, 
emitting or absorbing the well-known 21-cm(1.4 GHz) line due 
to the spin flip of the electron. This emission is redshifted by a 
factor 10-1000 due to the cosmological expansion and ends up 
in the frequency range from 140-1 .4 MHz. When the hydrogen 
spin temperature is not coupled perfectly to the radiation tem- 
perature of the cosmic background radiatior[3, but changed by 
other couplings with the surrounding matter and radiation, it 
can be seen against the cosmic background radiation in absorp- 
tion or emission, depending on whether the spin temperature is 
lower or higher than the background radiation temperature. In 
this way, the cosmological 21 -cm emission carries information 
about the evolution of the Universe. 

The exploitation of cosmological 21 -cm emission is sub- 
iect of a r apidly developing literature, which is reviewed by 



Furlanetto e t al. (2006). Here we highlight the requirements for 



of Reionization ( 93. 1.1b as well as from the Dark Ages be- 
yond reionization ( 93. 1.21 ). 21 -cm tomography of the reioniza- 
tion era, ( 93.1.31 1. and measuring the power spectrum of 21 -cm 
fluctuations ( 93.1.41 ) out to z=50. 



3.1.1. Global Epoch of Reionization 

The 21 -cm hyperfine transition of neutral hydrogen atoms 
can be used as a cosmological probe because differences be- 
tween the spin temperature characterizing the population of 
the upper and lower hyperfine state and the radiation temper- 
ature Tr of the background radiation are observable directly: 
for Ts < Tr, the 21 -cm emission is seen i n absorption, while 
it is s een in emission for the converse case dEwen and Purcell , 
195lj). The spin temperature is determined by the history of 
coUisional and radiative excitation, which in turn is determined 
by the cosm ological evolution of CMB and gas temperatures 



(see, e.g., [S unvaev and Z eldovichl 119721 : IScott and Reesl 



1990l) . lGnedin and Shaver (200^ have simulated the evolution 
of the sky-averaged brightness temperature of redshifted 21 -cm 
emission due to the onset of reionization, the so-called global 
EoR signal. The signal is a characteristic variation of 21 -cm 
brightness temperature T2\iz) with redshift, which arises at fol- 
lows. Before reionization, the hyperfine populations of neutral 
hydrogen are in equilibrium with the background radiation, 
and AT2\ = - Tj^(z) - 0. As the first Lyman-a photons are 
generated, becomes decoupled from the background radia- 
tion temperature and couples to the gas temperature, eith er b y 
collis i ons or via Wouthuysen-Field scattering dWouthuvsenl 



three apphcations: detecting the global signal from the Epoch 



Which is not a microwave background at those redshifts, of course. 



1952 : Field , 19591) . Since the gas temperature is less than the 
CMB temperature and therefore AT2\ < 0, the 21-cm line 
is seen in absorption. As the gas heats up, AT21 increases 
and reaches values greater than 0. The heating increases until 
reionization has proceeded so far that only very few neutral 
hydrogen atoms are left, thus removing the source of the 21- 
cm emission and resetting AT2\ to its initial value of 0. The 
observed wavelength at which this signal appears directly re- 
veals when reionization occurred, and the redshift evolution of 
Ar2i(z) encodes the detailed physics of reionization. 

Because of the homogeneity and isotropy of the Universe, 
this transition happened roughly simultaneously in all parts of 
the Universe. Therefore, the same signal is observed indepen- 
dently of the direction one is looking in, and a single simple 
dipole antenna would be sufficient equipment to pick up the sig- 
nal, which has an amplitude of a few to a few tens of milliKelvin 
dGnedin and Shaven , |2004 Fig. 1). Thus, although simple in 
principle, detection of the signal is difficult against the Galactic 
foreground. As reionization is believed to occur in the redshift 
interval 6 < z < 20, the signal will appea r at frequencies be- 
tween 70 and 200 MHz. As discussed by iGnedin and Shavei 
(l2004 . he main problem here is not the sensitivity as such. 
The sky temperatures in this frequency range are 100-2000 K, 
with lower temperatures at higher frequencies (eq.|5]). As Fig. [8] 
shows, the necessary observing times to reach 10-cr detections 
of a 1 mK signal from z = 15 with a single dipole are just a few 
months. Hence, the ability to distinguish the global EoR signal 
from the effect of the foreground variation imposes a further 
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Single-dipole global EoR exposure time 
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Fig. 8. Integration time necessary to reachanRMS noise level of 1 mK in the redshifted 21 -cm line for observing the global Epoch of Reionization signal (of 
order of a few to a few tens of mK: lGnedin an d Shaver, 2004, Fig. 1) with a single dipole against the Galactic foreground emission, as function of redshift 
(z = 5 corresponds to v = 230 MHz, z = 30 to 45 MHz) The integration time is obtained from eqns. |5] and with / = 1. A constant fractional bandwidth of 
5% was assumed. To reach fainter noise levels s, the required integration time scales as ?exp »: l/s^. The redshift range shown includes the entire likely range 
for the oiigin of the global EoR signal; redshifted 21-cm emission from higher redshifts is considered separately below. 



observability constraint. The evolution of Ar2i(z) corresponds 
to a smooth change of the spectral index as function of fre- 
quency, which needs to be distinguished from similarly smooth 
variations in the Galactic foreground emission. Its detection is 
- so far unsuccessfully - being tried from the ground and re- 
quires very high gain stability and absence of RFl: the Australia 
Telescope National Facility has an ongoing development ef- 
fort "Cosmological Reionization E xperiment" (CoRE) that at - 
tempts a globa l EoR measurement (iChippendale et all 120051) : 



Bowman et alj ( |2008|) have just reported first results from the 
EDGES experiments determining an upper limit of 250 mK 
to the reionization signal, and stressing the need for systematic 
errors being limited to less than 3 mK. 

A single antenna with an instantaneous frequency range 10- 
200 MHz, in an "eternally dark" lunar crater at the pole, kept 
at a stable temperature, and shielded from any solar and terres- 
trial radiation, would clearly offer the very best set-up for such 
an experiment. That would reduce most instrumental and RFl 
effects to the bare minimum and only leave the cosmic fore- 
ground as the main issue. Using several independent dipoles, 
without any beam-forming, would increase sensitivity and re- 
liability against systematic errors even further With enough 
signal paths the same antennas could then also be used in a 
beam-forming mode. Systematics and cosmic variance could 
be further checked by using both lunar poles which lets one ob- 



serve two independent patches of the sky. A lunar array could 
also collect the global EoR signal simultaneously with other 
observations by adding signals from different antennas incoher- 
ently. This only requires including a dedicated signal process- 
ing chain parallel to the correlator. Together this should deliver 
an ultra-precise global radio spectrum of the low-frequency sky 
to identify the brightness temperature change in the the 21 cm 
line caused by reionization. This would answer the question: 
When did the global transition between a predominantly neu- 
tral and a predominantly ionized Universe happen? 



3.1.2. The global 21-cm signal from the Dark Ages beyond 
reionization 

In the existing literature, the focus for 21-cm cosmology has 
been put on the global signal from the Epoch of Reionization on 
the one hand, and on the power spectrum of 2 1 -cm fluctuations 
from the EoR and higher redshifts on the other (see following 
sections). However, there is also a global signature of the dark 
ages from the spin temperature evolu tion at redshifts beyond 
reionization and up to z ~ 2 00 (e.g.. iCiardi and Madaul 12003 



" ihttp : //www, haystack .mi t ■ edu/ast/arrays/Edges/| 



Ciardi and Salvaterrai l2007h . The expected signal is a surface 
brightness decrement with a peak amplitude of around 40 mK 
at z a; 90 and approaching at z » 20 and z ~ 400. We have 
calculated the exposure time necessary to observe this global 
Dark-Ages signal with a si ngle dipole against the sky back- 
ground, using the model of Ciardi and Salvaterrai ( 2007 ). The 
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Fig. 9. A simulated global 21 - cm HI signal from the dark ages after 1 year of integration time with a single sky-noise limited dipole as a function of frequency 
(based on the calculations by ICiardi and Salvaterral , [20071 also compare Fig. llOt . Here a fixed bandwidth of 1.5 MHz is used. We assume that the observed 
brightness temperature, T^bs, is just the sum of the dark ages signal and a foreground signal, Tfai, which is a perfect power law. The signal at 30 MHz 
originates from z = 46. 




Fig. 10. Global 21-cm signal from HI in the dark ages. The plot shows the observing time needed to reach a 5-<t detection at 5% fractional bandwidth with 
a single dipole, as function of redshift. If N dipoles are added incoherently, then the observing time is reduced by a factor Vw. 
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result is shown in Fig. |9l with the exposure time necessary to 
reach a 5-cr detection for a signal peaking at different redshifts 
shown in Fig. [TO] A single dipole can provide a 5-cr detection 
at z = 42 in 137 days. A caveat here is that the presence of the 
reionizing Lyman-o' photons may be tightly correlated with high 
densities also beyond the very earliest stages of reionization. In 
this case, the final spin temperature may be much higher than 
the originally predicted diffuse hydrogen temperature, i.e., the 
spin temperature is closer to the CMB temperature, decreasing 
the absorption dip and weakening the absorption signature. 

The simulated line detection with Gaussian noise in Fig. |9] 
shows that the signal is only 10"^ of the foreground signal. 
Moreover, a detection has to rely on the foreground to be a 
perfect power law. Hence, the foreground subtraction is the 
biggest uncertainty together with the need for an almost per- 
fect band-pass calibration of the antenna. It is likely that the 
actual foreground signal will deviate from a simple power-law. 
This makes a detection almost impossible if there are spec- 
tral foreground fluctuations of the same magnitude as the dark 
ages signal. Nonetheless, further studies need to show on which 
level spectral foreground fluctuations are actually present. As 
in the global EoR case, the global dark-ages 21 -cm signal is 
isotropic, so that a single (well-calibrated) dipole suffices to 
carry out this measurement and detect the signal from Fig.|9]in 
frequency space. 



3.1.3. EoR tomography 

At the other extreme of the spectrum of potential cosmol- 
ogy experiments would be an attempt to image the hydrogen 
fluctuations directly. The detailed history of reionization can 
be traced by observing the two-dimensional structure of neu- 
tral and reionized gas around luminous objects, and by step- 
ping through the different frequencies corresponding to differ- 
ent emission redshifts, a tomographic map of reionization can 
be created. The characteristic scale of the fluctuations is ex- 



pected at the arcminute scale (Furl anetto et al.Ll2004^ . This adds 
arcminute-scale spatial resolution over the frequency range 50- 
140 MHz (corresponding to redshifts 10 < z < 20) to the re- 
quirements, as well as the ability to detect cosmological milli- 
Kelvin brightness fluctuations over the scale of one resolution 
element. The same hardware would allow to image hydrogen 
fluctuations directly also during the preceding dark ages, at 30 < 
z < 100, if the frequency range is extended down to 14 MHz 
and the baselines are increased correspondingly to achieve the 
desired resolution also at these longer wavelengths. 

The purity requirements for EoR likely mandate a lunar far- 
side location to ensure absence of solar, terrestrial auroral, and 
man-made interference. A near-polar region could be consid- 
ered since this will provide a near-constant view of the sky, 
allowing long integrations. If the site is closer than ~6.5° to the 
pole, lunar libration may expose the site to the Earth and make 
EoR observations impossible for some fraction of the year 

The expression for the exposure time necessary to reach a 
given surface brightness sensitivity Trms is particularly simple 
when expressed as function of the system temperature (which 
is a function of observing frequency, i.e., redshift, via V2i(z) - 



1.4GHz/[l -i-z]) and given by solving eq. (|7]) for the exposure 
time: 

^ ^ l_ /' rsky(v2l[z]) \' 



RMS 

2.1 X lO'^yr / b 
N{N - 1) \ 5% 



/TRMSrV^ 

llmK/ \1' 



-4 / , X 6.06 

1 + z\ 
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(25) 



where / is the usual filling factor (eq. [8]l and b - Av/v is the 
fractional bandwidth. Here, we have used the fact that the sky 
temperature happens to take on the simple form 



rsky(^2l[z]) = (l+z)"'xlK 



(26) 



using the approximation to rsky(v) from eq. (|5]l for the observ- 
ing frequencies v > 2 MHz considered here, and substituting 
V2\{z\ - 1.4GHz/(l + z) (see also upper panel of Fig. [T2]l. 

The exposure time from eq. $25[ is plotted against target 
redshift in Fig. [TT] assuming a scaled array in which the max- 
imum baseline increases with wavelength to keep the reso- 
lution constant at i? = 1'. This implies a filling factor / = 
2.1 X 10"^ -\/N(N - 1); conversely, the number of dipoles re- 
quired for a given filling factor at a given resolution is 



A^dipoles * 47 X 10" / l- 



-2 



(27) 



for A^dipoies s> 1 . The corresponding effective area is A^dipoies x 
^2/4, i.e.. 



A,o, = 10.5 km^ A^dipoies (—J Y7 



(28) 



where 1 MHz/v can be substituted for the /i/30m term. 

Again, the exposure time is driven by the large brightness 
of the Galactic synchrotron emission, which sets the system 
temperature. The exposure time can only be decreased by in- 
creasing the filling factor, i.e., the number of dipoles; since the 
filling factor is limited to / < 1, setting / = 1 in eq. ( l25t yields 
a firm lower limit to the necessary exposure time for a given 
Trms- In principle, if enough correlator capacity is available, 
the imaging could be performed over a significant fraction of 
the entire visible sky within that time scale, since dipole anten- 
nas have a very large instantaneous field-of-view (compare the 
discussion of survey speed in ^2.5l l. 

It is evident that an array with a substantial number of dipoles 
is required to reach the sensitivity necessary for 21 -cm tomo- 
graphic imaging. For observing times of order one year for 
imaging around z = 20 (at 66 MHz or 4.4 m), of order 10^ 
dipoles would be required, with a total effective area of 50 km^ 
(also see Fig. |4|. The strongest factor here is the spatial res- 
olution which enters with the fourth power in the integration 
time (eq. \25[ and squared in the number of dipoles needed 
(eq. l27i. Hence, a 10' a rray observing 21 -cm emission from 
z = 20 (IZahn et al.Ll2007h would require only about 10^ dipoles 
(i.e. 0.5 km^) at the same sensitivity limit. Higher redshifts - 
in the pre-EoR area - become increasingly more difficult even 
though the collecting area per dipole increases with (i.e., 
(1 -H z)^) because the sky temperature increases rapidly with 
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Array for high-z 21cm tomography 




z 



Fig. 11. Number of crossed dipoles (from eq. (3]) required for a 3-o" detection of 1 mK per resolution element in one year when imaging the redshifted 21cm 
emission at frequency V2i(z) = 1.4GHz/(l +z) at angular resolution & = 1', assuming a constant filling factor independent of wavelength (i.e.. a scaled array). 
The required number of crossed dipoles to achieve this RMS can be traded against longer exposure times (and vice versa) as indicated, but must not exceed 
the shadowing limit from eq. ^9)' ^ < ^max = 4.7 X 10^(??/1') "■ The con'esponding collecting area is given by 1.1 X 10"** ^(1 + z)^ km^. The an'ay diameter 
(or maximum baseline) required to reach a given resolution is shown in Fig|5]and ranges from a: 800m (z = 10, if = 10') to » SO km (z = 100,!? = 1'). The 
numbers from this graph can be scaled as Nam V'^p ' "^"(^ = const (where S is the desired sensitivity). 



longer wavelengths. Fig. Hand eq. dZSl l imply that at the highest 
redshifts the number of dipoles needed would require a filling 
factor larger than unity for one year of integration time and 1 ' 
resolution. This is of course impossible and can only be reme- 
died by using antennas with a higher gain than the crossed 
dipoles assumed here. For example, log-periodic or Yagi an- 
tennas have a higher gain (i.e., eff'ective area) per antenna, at 
the cost of a smaller field-of-view and, worse, increased com- 
plexity compared to simple crossed dipoles. As Fig. [TT] shows, 
a filling factor / = 1 is already needed to achieve exposure 
times of order 1 year for tomography at z x 15; since the ex- 
posure time at fixed resolution and filling factor grows roughly 
as (1 -H z)^, substantially higher redshifts can only be reached 
with very high-gain antennas or very long observing times. 

Confusion noise poses further constraints: 21 -cm emission 
from z = 15 is observed at 87.5 MHz, where the confusion 
noise at 1' resolution is about 14mJy, which corresponds to a 
brightness temperature of 710 K. This is much more than the 
desired signal of 1 mK. To resolve the confusing background 
and reduce it to the 1 mK level, which is prerequisite for allow- 
ing its subtraction, it would be necessary to improve the reso- 
lution to the (X'Ol level, which would require baselines larger 
than the diameter of the moon (3476 km). Instead, again statis- 
tical techniques have to be employed to extract the signal from 
the confusion noise via their different structures in frequency 
space (iDi Matteo et al..i2004) . 



In summary, one can conclude that a lunar EoR imaging ex- 
periment needs to be very substantial in size to penetrate into 
the z > 20 epoch, but is not unthinkable in the long run. This 
crucially depends in the expected angular size scale of the sig- 
nal, which will hopefully become more clear with advances in 
ground based experiments. Accurate subtraction of the sources 
of extragalactic confusion noise is also a necessity and may 
prove difficult. 



3 . 1 .4. Power sp ectrum of the 21 cm-transition at z — 30-50 

Loeb and Za ldarriaga (2004) have pointed out that the fluc- 
tuation power spectrum of the redshifted 21 -cm line carries a 
wealth of information about the matter power spectrum. The 
CMB radiation itself that is observed today carries information 
about cosmological parameters mainly at the largest angular 
scales from 0?2-90° (multipole scale I from 2-1000, and lim- 
ited in principle by Silk damping to I < 3000). By contrast, 
the angular power spectrum in the redshifted 21 -cm line car- 
ries cosmological information at much smaller angular scales, 
I > 10"*, corresponding to angular scales of V or less. In addi- 
tion, observations of redshifted 21 -cm emission from different 
redshifts in the range 30-50 yield independent samples of the 
cosmological parameters, while the CMB information suffers 
from cosmic variance. Therefore, observations of the 21 -cm 
power spectrum from this redshift interval yield many orders of 
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magnitude more information about density fluctuations in the 
early Universe than direct observations of the CMB power spec- 
trum. This makes the redshifted 21 -cm line a powerful tool to 
constrain all model parameters necessary to describe the Uni- 
verse, e.g., the slope and curvature of the initial density fluc- 
tuation spectrum containing information from the inflationary 
phase, the mass of warm dark matter particles, and the fraction 
of matter density contributed by neutrinos. 

In particular, angular scales I > 5000 carry the greatest 
weight in discriminating between different values for the cos- 
mological parameters, since it is at these wavenumbers where 
the predicted power sp ectra for d ifferent sets of cosmological 
parameters differ most CLoeb and Z aldarriaga, 2004, Fig. 3). A 
low-frequency measurement of the 21 -cm power spectrum at 
these angular scales and in this redshift interval would there- 
fore be an ultimate cosmological experiment. 

Baryon acoustic oscillations (BAOs) are a major tool of 
modern cosmology, producing a feature that is directly vis- 
ible in an observed power spectrum. They are imprinted on 
the 21 -cm angular power sp e ctrum in the same way as on 
flie CMB (de Bernardis et al. ', '2000^ and galaxy po wer spec- 
tra fe isenstein et al., 2005; Cole et al., 2005; Perci val et all 
l2007h . The BAO sig nal is present on spatial s cales of 20- 



400 Mpc (comoving; ' Barkana and Loebl 120051 Fig. 4). In 
the 21 -cm angular power spectrum, it is predicted to appear 
at wavenumbers of order / « 90^000 and with oscilla- 
tion a mplitudes of order or milliKelvin or smaller (se e, e.g., 
Lewis and ChalHnor, 2007, Fig. 8 and iMao and Wul |2008, 



Fig. 4). 

The distinct advantage of observing the 21 -cm power spec- 
trum from the Dark Ages, before the Epoch of Reionization, 
is that the astrophysics of reionization is complex, non-linear 
and not well understood, but these complexities do not affect 
the Dark- Ages power spectra, which therefore provide a much 
cleaner signal to observe and to interpret. For exa mple, the non 



Gaussianity of the fluctuations can be probed (ICooray et al. 
[2008) and may reveal information on the inflationary phase of 
the universe. 

The required observing frequency is 1.4 GHz/(l + z); see 
the upper panel of Figure [121 As noted in section 12.4.51 angu- 
lar scatter broadening in the ISM and 1PM limits the highest 
observable angular scale as a function of frequency. The lim- 
iting angular scale increases proportionally to (see Fig. |5] 
for the maximum possible resolution and the baseline neces- 
sary for reaching this resolution). The lower panel of Figure[T2l 
shows the number of independent samples (data points) that 
can be collected up to the IPM scattering limit (A^2icm from 



Loeb and ZaldarriagaLl2004 p. 4), assuming a bandwidth of 5% 
of the observing frequency. Since the number of samples scales 
as /^j,^ and /max scales as A~^, which in turn scales as the 
number of samples up to the scattering limit scales as (1 +zy^- 
Therefore, observing fluctuations originating from z = 45 only 
contributes 10% of the information accessible at z=30, and ob- 
serving z = 70 only contributes an additional 1 % compared to 
z = 30. Thus, as far as the number of samples is concerned, 
the presence of angular broadening in the ISM and 1PM makes 
extending this experiment to lower frequencies progressively 



more inefficient. However, the detailed strength of the signal as 
function of redshift cannot be predicted in advance, and there 
may be highly relevant deviations from current theory. This 
is particularly important since, as Loeb and Zaldarriaga (2004|) 
point out, this phase of the cosmological evolution is a fairly 
clean laboratory, hardly affected by any 'dirty' astrophysics of 
the reionization process. Measuring the evolution of structure 
formation with redshift between z = 10 - 100 may therefore 
be one of the most intriguing diagnostics of the history of the 
Universe. For these reasons, it will still be desirable to cover a 
broad redshift, i.e. frequency, range. 

Sensitivity poses another constraint. The angular power 
spectrum has the strongest total signal at z ~ 55 (25 MHz), 
dropping off roughly symme trically to z = 25 and z = 170 



(ILoeb and Zaldarriaga! 12004 Fig. 2, e.g.). The signal strength 
is expected to be at the mK scale, while the system temperature 
and hence the noise are determined by the sky temperature (see 
Figure O. Figure [13] shows the observing time necessary to 
reach an RMS level of 1 mK in the power spectrum wi th an ar- 
ray of dipoles as given bv lLoeb and Zaldarriagal lEooi eq. 11). 
To achieve an angular resolution of 1' in reasonable integration 
times of order one year, about 10^^ dipoles are needed for z = 
30 (Aeff =3.5 km^) and 10^ dipoles (Aefr =30km2) at z = 50. 
Again, this number depends very critically on the angular res- 
olution required and would reduce to N ~ 10^^ and * 10"^, 
respectively, for angular scales of 10' if the signal were on the 
mK level. Distinguishing between different cosmological mod- 
els at high confidence levels or detecting BAOs directly will 
probably require sensitivities down to the 0.01 mK level, but 
the latter may already be feasible with degree-scale resolution, 
probing multipole numbers I x 100 (Barkana and Loeb, 2005). 

Again, accurate subtraction of spectral and noise power 
fore- and backgrounds is essential. Determining the exact ar- 
ray parameters that will be necessary requires more detailed 
guidance from further cosmological observations and simula- 
tions. Numerous novel approaches to extracting cosmological 
information from 21 -cm power spectra are currently being 



iMorales et al., 20061; 


Barkana and Loebl 2008 ; 


Visbal et al. 


,2008; Cooravetal., 


|2008|). 





3.2. Extragalactic surveys: radio galaxies, clusters and 
feedback 

At low radio frequencies, the brightest sources of the extra- 
galactic sky are the synchrotron-emitting radio lobes of giant 
radio galaxies and quasars, such as Cygnus A. In fact, the 178 
MHz 3C survey ( now used in its twice-revised form 3CRR, 



Laing et all Il983h was the first successful survey of discrete 
sources, and 3C sources are still studied extensively today. 
There are a number of questions that are currently of particu- 
lar relevance to which the study of radio galaxies can provide 
decisive insights. 
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Fig. 12. Upper panel: observed frequency of redshifted 21cm emission (left-hand axis), and baseline required to reach 2' resolution at that frequency (right-hand 
axis), both as function of redshift z- Ground-base d observations are exce edingly difficult below 30 MHz, corresponding to z = 47. Lower panel: Number of 
independent power-spectrum samples (A'2icm from iLoeb and Zaldarriagal . l2004l p. 4) obtainable down to the ISM scattering limit towards the Galactic poles, 
also as function of redshift. 
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Fig. 13. Number of antennas necessary to achieve a 5-cr detection of fluctuations in the 21cm power spectrum in one year at multipole numbers up to 
/ = 50 00 (corresponding to 2' angular resolution) or / = 500 (20') at the 1 mK level, as function of redshift (as given by eq. 11 from Loeb and Zalda rriag^ , 
l2004h . Except for the dift'erent exposure time prescription, this figure is constructed in the same way as Fig. so that the same scalings apply, i.e., 
^Ani V^e-^p ' '^(1 = const (where S is the desired sensitivity level), assuming a filling factor independent of frequency. 
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3.2.1. Electron energy distributions and jet power in radio 
galaxies 

The low-frequency emission of radio galaxies and quasars 
is generated predominantly in the "lobes". These are struc- 



tures with scales of t ens to hundreds of kpc dMilevl \l9?,(i 
iBegelman et al. , 1984 ) that are inflated by pairs of relativistic 
jets emerging from near the accretion disk around a black hole 
in the nucleus of the host galaxy. 

Jets and lobes of radio galaxies emit synchrotron radiation, 
which has the property that the lowest-frequency emission is 
generated by the lowest-energy electrons. Since electron energy 
distributions typically are of the form N{E) dE ~ E^p AE with 
p a: 2.4 to yield the observed spectra v"*^^, low-energy electrons 
are the most numerous and dominate the power output. On the 
other hand, as discussed below, higher-energy electrons have 
shorter cooling time scales, leading to a steepening of the radio 
spectrum - the so-called spectral agi ng ([Alexander and Leahvi 
1987 : Blundell and Rawlingsl 200lV Hence, older parts of a 



source tend to have steeper radio spectra. Low-frequency ob- 
servations thus probe the oldest structures in a galaxy and can 
be used to constrain the age of a source. 

If a source is sufficiently compact, it can become optically 
thick to its own synchrotron emission due to synchrotron self- 
absorption. This changes the observed spectral energy distri- 
bution from the intrinsic optically thin synchrotron spectrum 
fv oc v"", with a X -0.7 as typical value, to fy oc y^^^, censor- 
ing some of the observables by hiding spectral structures and 
effectively removing access to the low-energy electrons. The 
self-absorption limited maximum brightness t emperature scales 
asjrnax ~ lO'^ K (v/ 1 MHz)!/^ (fi/ inT)- (IBegelman etal. . 
19841 eq. A8). Typical magnetic fields in lobes are estimated to 



be of the order of a few nanoTesla. Lobe surface brightnesses 
are wen below lO'^K, the maximu m allowed by "inve rse- 
Compton cooling catastrophe" ( Kelle rmann and Paulinv-T oth, 
The lobes of radio galaxies are typically not affected by 
optical-depth effects. 

This is, however, not true for the "hotspots", the terminal 
shock front of powerful jets. These frequently show a spectral 
turnover at frequencies below 100 MHz. From the size of a 
hotspot and the turnover frequency one can deduce reasonable 
limits of the magnetic field strength and power of the jet feeding 
it. 

The inflation of lobes by jets means that mechanical work 
is done on the surroundings of the radio galaxy. Jets and lobes 
therefore provide a means for a mechanical coupling between 
accreting black hole at the centre of the radio galaxy and its 
surroundings. In current mo dels of galaxy and black hole for 
mation and evolution (e.g. iBower et al. , l2006t ICroton et al 



2006t iHopkins et all l2007h . this AGN feedback is key to 



regulating the growth of the black hole and establishing the 
locally observed correlation between black-hole mass and 
large-scal e galaxy properties such as bulge mass and bulge lu- 
minosity ( Ferrares emd Merritt , I2OOOI: iGebhardt et al.[ . I2OOOI: 
Haring and Rix. 200^ . 



To calibrate and verify these galaxy formation models and 
to understand the ICM thermodynamics, it is very important to 



determine the magnitude of this mechanical energy input into 
a radio galaxy's surroundings. Since the total energy of the 
radio lobes is dominated by the lowest-energy electrons, low- 
frequency observations are essential. Yet, today, still very little 
is known about the shape of the electron energy distribution at 
frequencies below 10-100 MHz. 



3.2.2. High-redshift galaxies and quasars 

One successful method to find high-redshift radio galaxies 
quasars i n the past was to search for ultra-steep-spectrum radio 
sources (^ Roettgering et allll997 ). Since observations at lower 
frequencies include preferentially steeper-spectrum sources, a 
ULW survey might turn up a large number of high-redshift ra- 
dio sources. Extrapolating (perhaps optimistically) from source 
counts at 74 MHz in the VLSS (eq.l2TI). our example array from 
the same section (300 elements, 100 km maximum baseline, 
10 MHz observing frequency) would thus be able to discover 
3 Million sources down to the confusion limit of 65 mJy at 5 o" 
significance over half the sky (500000 per steradian), which it 
would be able to survey within roughly two years. For refer- 
ence, Roettgering et al. ( 1997) found 30 radio galaxies at z > 2 
starting from a sample of ^ 700 ultra-steep spectrum sources. 

Another, not yet well-explored method to identify active 
black holes and jets in the early universe could be to look 
for young jets that are still contained within the just forming 
galaxy. They are expected to show up as compact radio sources 
with synchrotron self-absorption turnover frequencies that are 
unusually low for their sizes (lFalckeetaUl2004 . In the local 
universe these sources are known as Compact-Steep-Spectrum 
(CSS ) and GHz-Peaked-Spectrum (GPS) sources (O'DeiJ 
with turnover frequencies in the range 100-1000 MHz. 



At z = 10 these would be shifted to 10-100 MHz and could 
be readily identified by an ULW array. Interestingly, the only 
radio-l oud quasar so far identified at z > 6 may be a CSS 



source (IFrey et al.Ll2008h . 



The number-counts estimate above, of course, includes all 
kinds of radio sources, so that the 500000 sources per stera- 
dian will also include objects from this category and the ones 
described next. 



3.2.3. Fossil radio galaxies 

Radio galaxies and other AGN are now believed to be nor- 
mal galaxies in a phase of more or le s s stron g accretion onto 

20061) . It is an impor- 



their central black hole ( ICroton et al 



tant question how long this active phase lasts, and on which 
timescale their activity recurs. 

In the so-called "double-do uble" radio galaxies with tw o or 
even sets of radio lobes (e.g. Schoenmakers et al. , 2000l) . the 
recurrence is directly visible. From synchrotron ageing argu- 
ments, the age of the outer set of lobes is inferred to be of 
the order of several tens of millions of years, while the inner 
lobes are inferred to have been launched after an interruption 
of several to several tens of millions of years. Thus, both the 
duration of an activity phase and the recurrence timescale are 
likely within an order of magnitude of 10^ yr 
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However, double-double sources are rather rare, perhaps im- 
plying that special circumstances might be required for their 
formation. Other methods are needed to determine the duration 
of a galaxy's active phase in the general case. A very direct 
way of doing so is counting how many "fossil" radio galaxies 
there are, i.e., those galaxies whose activity phase has switched 
off recently. However, the fossil galaxies have been elusive 
so far, in part because the emitted power of a radio galaxy 
drops as it grows in linear extent even for constant jet power 
dBlundell et al., 1999). This creates a "youth-redshift degener- 
acy" in current flux-limited low-frequency radio surveys: only 
fairly young radio galaxies (ages ~ 10^-10^ yr) are observ- 
able at high redshifts (z > 0.8), while at lower redshifts, radio 
sources are observabl e in principle up to ages of 10^-10"^ yr 
( Blundell et al. . 199 9^. Fig. 17). Because of this degeneracy no 
strong statements can be made yet about radio source lifetimes 
based on the study of low-frequency samples. 

Moreover, the only known radio galaxy relicts have been 
identifie d as cool holes in the X-ray emission o f the h ot cluster 
gas (e.g. iMcNamara et all I2OOII: iFabian et all l2006i) . not via 
direct searches for their radio emission (see also the following 
section on galaxy clusters). Present observations simply may 
have been done at frequencies that are too high to detect fossil 
radio galaxies in high-frequency all-sky surveys. This is likely 
to be the case also with the much better sensitivities of the up- 
graded E-VLA and E-MERLIN arrays. Lower-frequency ob- 
servations are better suited for detecting relict radio sources di- 
rectly because they probe electrons with lower energies, which 
are those with the longest radiative loss times. This is because 
of two key properties of synchrotron radiation: 

(i) The frequency at which a relativistic electron of energy 
£■ in a magnetic field of flux density B emits synchrotron 
radiation scales as 

Vc oc E^B, 

i.e., lower-frequency emission is generated by lower- 
energy electrons. 

(ii) The energy loss of a synchrotron-emitting electron scales 
as 



AE , , 

oc E^B^, 

At 



(29) 



i.e., synchrotron cooling times are inversely proportional 
to energy. 

But it is not only the case that synchrotron cooling times in- 
crease with decreasing observing frequency - more strongly 
yet, there is a maximum energy that a synchrotron-cooling elec- 
tron can have after a given time has elapsed since its accelera- 
tion e ven if it had been given an infinite amount of energy ini- 



2001), be 



tially (Ivan der Laan and Perolal Il969t iJester et al 
cause eq. (|29T l integrates to 1/E{t) - l/E{to) oc f - fg. 

Hence, observations at lower frequencies than currently ac- 
cessible allow fossil radio galaxies to be detected for a longer 
interval between the switching-off" time and the time of ob- 
servation, with the accessible radio-lobe lifetime increasing 
roughly inversely proportional to the observing frequency. Low- 
frequency radio observations may even be the only possibil- 



ity for detecting relict radio galaxies. Even the continued non- 
detection of such fossil galaxies would yield new physical in- 
sights, since it would put the strongest constraints on how 
rapidly radio lobes have to dissipate all of their energy after the 
central engine ceases to be active. 

3.2.4. Galaxy clusters 

Just like the gas in radio galaxies, the plasma that fills galaxy 
clusters emits synchrotron radiation, as well as X-rays through 
inverse-Compton scattering. The origin of the radiating rela- 
tivistic particles is unknown. Theories for their origin can be 
constrained by low-frequency observations, analogous to those 
of radio galaxies. This will yield important constraints on the 
formation history of the clu sters and of their magnetic fields. 
Cassano et al.l ( 20061 12OO8 ) have shown that the steep spectra 
of cluster halos mean that observations at lower frequencies 
are much more likely to detect the cluster hal os, specifically 



for lo w-mass and high-redshift objects. Indeed. iBrunetti et al 



(12008 ) report the detection of diffuse steep-spectrum emission 
from the cluster Abell 521 at 240-6 10 MHz that is not de- 
tectable at 1 .4 GHz. Again, we stress that even a fairly moderate 
array of 300 dipoles and 100 km maximum baseline operating 
at 10 MHz would be able to discover hundreds of thousands of 
sources in a matter of a few months. 

Independently of the question of the impact of AGN feed- 
back on galaxy evolution discussed in ^3.2.11 above, the me- 
chanical energy input into the intra-cluster medium (ICM) from 
giant radio sources is an important quantity in studying the 
thermodynamics of galaxy clusters, and in erudating how en- 
ergy is transported through them in order to establish the ob- 
served temperature and entropy profiles. For example, the de- 
tailed X-ray observations of the Centaurus and Perseus clus- 
ters bv iFabian et al.l ( 120051 120061) show that bubbles blown by 
the central radio source play an important role. In the case 
of the Perseus cluster, the X-ray images have revealed sets of 
"ghost bubbles" without detectable high-frequency radio emis- 
sion, which are identified as relict sources as discussed in the 
previous section. These ghost bubbles rise buoyantly in the 
cluster atmosphere and are therefore part of the energy trans- 
port mechanism in clusters. With low-frequency observations, 
it may be possible to detect radio emission from known bub- 
bles, or new bubbles in clusters that have not or cannot be 
observed to similar X-ray depth as Centaurus and Perseus. In 
either case, low-frequency radio observations would yield the 
most stringent constraints on the ages of such bubbles. 

3.3. Galactic/ISM surveys 

In this section we discuss how to infer the structure of the 
interstellar medium and the origin of cosmic rays through low- 
frequency observations. Beyond these, also the study of pulsars 
may benefit from these because a confusion limit of 65mJy 
at 10 MHz with 1' resolution would be matched to the flux 
limit of LOFAR for the v"^ spectrum typical of pulsars. This 
implies that pulsars could either be detected directly, or limits 
can be set on their low-frequency cutoff. However, the impact of 
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temporal broadening in the ISM 02.4.6I I on the low-frequency 
observability of pulsars needs to be assessed in more detail than 
we can achieve here. 

3.3.1. Structure of the ISM - the Solar System's 
neighbourhood 

The structure of the ISM in general, and of the supernova- 
inflated "local bubble" within which the solar system resides, 
has been inferred by considering absorption li nes in u ltravi- 
olet and X-ray spectra of stars and the ISM (F^ischl 



Breitschwerdtll200UlOegerle et al.L l2005UHenley et al 



1998 



2007 



e.g.) and t he distributio n of p ulsar dispersion measures in the 
sky ( Cordes and Laziol 2002h . 



Also low-frequenc y radio observations have been brought to 
bear on this problem: P eterson and Webben (120021) have mod- 
eled the dumpiness of the free electron distribution in the warm 
ionized medium (WIM) phase of the ISM by matching the ex- 
pected turnover in the spectrum of the Galactic synchrotron 
e mission t o that observed by the RAE-2 (see ® and IMP- 
6 (IBrownl Il973h satellites. Their method models the radiative 
transfer through the WIM and requires only knowledge (or rea- 
sonable assumptions) about the intrinsic, unabsorbed shape of 
the underlying synchrotron emission spectrum in addition to 
the observed brightness distribution. It is applicable over the 
frequency range where the WIM is neither completely transpar- 
ent nor completely opaque, i .e., from 0.1-lOMHz according to 
Peterson and Webbed (l2002h . 

This modeling can be extended to infer not only the dumpi- 
ness of the ISM, but also the full three-dimensional structure 
of the WIM. Figure [14] shows the dependence of the ISM "vis- 
ibility" (distance to the t - \ surface) as function of Galactic 
latitude for different frequencies, and as function of frequency 
for different wavelengths. The limiting power-law for the line- 
of-sight distance as function of frequency (right-hand panel in 
Fig. [T4l i is given by the inverse of the free-free-optical depth 
(eq. fTTT i for the electron density and temperature in the mid- 
plane of the Galaxy; for the parameters given in the figure and 
at wavelengths between 0. 1 and 1 MHz, we obtain 



^^^=^^^^00^^013^) (tMH^)-^^«^ 
By scanning through these frequencies and modelling the ob- 
served emissivity as function of Galactic coordinates, the struc- 
ture of the WIM in the solar system's neighbourhood can be de- 
termined in 3D. The requirements on resolution and sensitivity 
are rather modest: taking lOOpc as a typical ISM bubble size 
and 1 kpc as a typical distance leads to a degree-scale resolu- 
tion requirement. The sensitivity requirement is comparatively 
modest also, since the ISM emission itself sets the background 
temperature, so that a sensitivity of order of one part in 1000 
of the background level is sufficient, i.e., the requirement is 
Trms ~ lO'^Knear 1 MHz. 

Around the turnover frequency of 2 MHz, the sky is uni- 
formly foggy; below it, the t - \ distance scales as given by 
eq. (l30b . Above it, the ISM becomes transparent at lower and 
lower angles above the midplane as the frequency is increased, 
i.e., the turnover frequency decreases with the angle above the 



>3/2 



midplane. Around 2.5 MHz, the entire Galaxy is transparent, 
while at 0.3 MHz, a frequency well above the putative lunar 
ionosphere, one would see out only to a few tens of parsec, and 
see the Galactic plane in absorption against the brighter WIM 
emission. At these low frequencies, the visibility will therefore 
be set by the actual substructure of the ISM (supernova shells 
and bubbles) rather than by the global structure. 

On the upside, the short distances out to which we can see at 
low frequencies imply the effects of ISM angular and temporal 
broadening are much less severe than for extragalactic obser- 
vations a^ii^he^freguencies^or example, from the NE2001 
model ( Cordes and Laziol 2002), the pulse broadening time at 
1 MHz over a distance of 500 pc in the direction of the galactic 
pole is about 90 s, instead of a few days for extragalactic ob- 
servations at the same frequency. Similarly, the angular scatter 
broadening is 0.5' instead of » 30', affording a linear resolu- 
tion of 0.1 pc at the t = 1 distance of 500 pc (see Fig. [TSlbdow 
for an indication of how angular and temporal broadening scale 
with distance and frequency). 



3.3.2. Origin of cosmic rays 

Following the overview by iDuric (l2000l) . this section de- 
scribes an approach for elucidating the origin of cosmic rays, 
i.e., relativistic electrons and protons, with the help of a low- 
frequency radio array. 

The essence of the approach suggested by 'Durid is to make 
use of the fact that Galactic HII regions are optically thick at 
frequencies below about 30 MHz. Therefore, radio emission 
observed towards such HII regions at low frequencies predomi- 
nantly arises from material along the line of sight to the HII re- 
gion, i.e., the synchrotron emission from cosmic-ray electrons. 
By observing the synchrotron emissivity towards different HII 
regions at different distances from the Earth, a 3D map of the 
electron density distribution can be built up. To remove the de- 
generacy in the synchrotron emissivity between electron den- 
sity and magnetic field strength, the radio data can be compared 
with gamma-ray observations. 

Supernova remnants (SNRs) are known to accelerate high- 
energy particles through the first-order Fermi shock accelera- 
tion mechanism. However, the so-called "injection problem" 
has not been solved, as the Fermi mechanism can boost parti- 
cles which are akeady mildly relativistic to very high energies, 
but it cannot accelerate thermal particles to relativistic energies. 
Since the frequency of synchrotron emission scales with the 
particle energy, low-frequency radio observations can be used 
to trace the energy distribution of the lowest-energy particles. 
This is necessary to constrain theories which attempt to solve 
the injection problem. 

The level of turbulence generated by SNRs in the interstel- 
lar medium is also relevant to the generation and propagation 
of cosmic rays. By observing the scatter broadening of extra- 
galactic background sources at various distances from a SNR, 
the level of turbulence can be probed and mapped in 3D. 

However, the resolution required to clearly resolve HII re- 
gions and SNRs is of order 1", and the required sensitivity 
is 0. 1 mJy/square arcsecond. Both are a challenge for a lunar 
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Fig. 14. Left: Distance to r = 1 surface of the Wann Ionized Medium component of tlie ISM due to free-free absorption, shown as function of Galactic latitude 
for different frequencies. The figure assumes the exponential electron density profile shown and neglects any azimuthal structure of the ISM. The radius of 
the circle is lOkpc, and the galactic plane is along the horizontal. The density at the plane was chosen to match the transition betwee n total absorption and 
unabsorbed sight lines in the polar directions to the observed turnover frequency of the ISM emission; for a more detailed approach, see lPeterson and Webbed 
(2002). The plotted profiles will be modified in practice by a higher electron density towards the Galactic centre, and by taking into account the detailed 
substructure of the ISM. Right: Distance to r = 1 surface assuming the same azimuthally symmetric exponential electron density profile as in the left-hand 
panel, now as function of frequency for dift'erent angles above the midplane. Around the turnover frequency of 2 MHz, the sky is uniformly foggy; at lower 
frequencies, the galactic plane is seen in absorption against the brighter high-latitude emission, while above the cutoff, the ISM becomes less foggy, beginning 
at high frequencies and high galactic latitudes. The limiting power law is given by D(t = 1) k 500pc(v/MHz)2 (see eq.|30). A gain, substructure will modify 
the r = 1 distance along different lines of sight. 



array and may be better done by a ground-based telescope at 
higher frequencies and resolutions. 

3.4. Transient radio sources: suns and planets 

Transient sources are particularly interesting in astrophysics 
because they often indicate very high-energy phenomena, e.g.. 
Gamma-ray bursts (GRBs), supernovae and other stellar ex- 
plosions, outbursts of accreting black holes, neutrons stars and 
white dwarfs, reconnection events in the solar corona, or the im- 
pact of electrons from the solar wind on planetary atmospheres. 
Other sources show transient activity simply because of rota- 
tion, most notably radio pulsars. Finally, at low frequencies the 
possibility for detecting coherent emission processes increases 
due to the macroscopic wavelengths. To allow the detection of 
transients, high-resolution high-sensitivity observations have to 
be possible within the timescale of the transient phenomenon. 

3.4.1. Solar and planetary radio bursts 

Both the Sun and planets produce radio bursts, the sun 
through magnetic activity in its corona, and the planets and 
their moons through interaction of charged particles in the solar 
wind with the planet's or moon's magnetic fields, as is the case 
for exoplanets (see section I3.4.2l i. In the case of the Sun, the 
most violent coronal events give rise to coronal mass ejections 



(CMEs), disruptions in the solar wind in which large numbers 
of charged particles (i.e., cosmic rays) are accelerated. The 
magnetic activity and the charged particles can have dramatic 
effects on the Earth and Moon, posing a danger to astronauts 
and electronic equipment outside the Earth's atmosphere, or 
disrupting radio communications. These effects have been 
called "magnetic storms", and efforts to predict the strength 
and timing of CMEs impacting on the Earth-Moon system are 
known as "space weather prediction". The CMEs themselves 
are traceable via their low-frequency radio emission as they 
propagate towards the Earth with an array capable of per- 
forming imaging at high time resolution. Since the emission 
frequency is dependent on density and thus radius, observa- 
tions down to 100 kHz would trac e emission out to ab out half 
an AU (using the density model of Mann et al. . 1999[ e.g.). 

While the Sun is the strongest source of radio emission at 
high frequencies (at 100 MHz and above), bursts from the mag- 
netized planets, in particular Jupiter, can reach a strengt h com- 
parab le to that of the solar bursts or more, i.e., a few MJy dZarka , 



2007t). The emission region of Ju piter is very small and not re 



solved on lunar-like baselines (see lNigl et al.L 120071) . Hence, an 
imaging application for Jupiter radio emission is not a scien- 
tific driver for a lunar telescope. However, since the emission is 
strong enough, a single broad-band receiver covering frequen- 
cies < 30 MHz for spectral measurements with high frequency 
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and temporal resolution of solar and Jupiter bursts could be 
envisaged for prototyping and background measurements. For 
angular discrimination between different burst sources, a small 
interferometer with comparatively modest baselines (beginning 
at 500 m for the highest frequency of 30 MHz) would be suf- 
ficient. In the long term, also interferometry on a Moon-Earth 
baseline could be considered to constrain Jupiter bursts in new 
ways. 



3.4.2. Extrasolar planets 

While Jupiter's emission is an additional foreground and 
interference signal that likely needs to be avoided for high- 
sensitivity low-frequency observations, it implies that low- 
frequency radio observations offer the most favourable achiev- 
able intensity contrast between an extrasolar planet and its sun, 
and therefore open up the possibility of a direct detection of 
extrasolar planets and their magnetosphere. 

IZark a (2007) has reviewed the radio emission from extra- 
solar planets and the results of searches to date. Magnetized 
planets emit bursts of low-frequency radio waves through the 
cyclotron maser instability (CMI), in which electrons from the 
stellar wind interact with the planet's magnetic field. The CMI is 
100% circularly polarized, while solar/stellar bursts are usually 
unpolarized, offering a mechanism of distinguishing between 
stellar and planetary radiation. Without a polarization-sensitive 
array, planetary and solar bursts would be distinguishable by 
temporal structures, either within burst s or from modulations 
on the timescale of the planetary orbit. ILazio et al](l2004l) give 
15 minutes as typical burst duration, while I Zarkal ( 200 7^hows 
an example of a burst with an overall duration of several hours, 
with frequency structure down to the millisecond level (in his 
nomenclature, only the features on sub-second timescales are 
called "bursts"). 

Since the CMI is only active in the presence of strong mag- 
netic fields, either in the planet or the stellar wind, this mech- 
anism is most suitable for the detection of magnetized Jupiter- 
like planets. Earth-like planets will only give rise to radio bursts 
if they are embedded in a stellar wind which is itself strongly 
magnetized. Even if this is a rare situation, the radio emission 
might be the most readily detectable signal of such planets. 

A crucial parameter is of course the rate at which plan- 
etary bursts occur, and hence the survey time that needs to 
be invested in order to obtain a statistically relevant sample 
of burst observations. However, very little is known about 
the likely recurrence rate, which depends on the unknown 
way in which CMI bursts are triggered. One possibility is that 
they have "light house beams" similar to pulsars, with intrin- 
sically continuous emission that is confined to a solid angle 
and bursts being obser ved when the emis sion beam sweeps 



over the telescope (e.g., Za rka et al.L 120041) . Another possibil 



ity is that CMIs are trigge red or enhanc ed by coronal mass 
ejections (CMEs; see lOurnett et al. . 2002 , who give an exam- 
ple of a solar CME enhancing Jovian radio emission). In this 
case, the relevant timescale is the CME recurrence timescale 
modified by the probability that a CME i n fact interacts with 
the planet. Based on such considerations, iKhodachenko et al 



(|2007b expect CME rates in excess of the solar value of about 
6 CMEs/day, in particular for M dwarfs which are intrinsic ally 
more active than solar- type stars. ' GrieBmeier et al. (2007a ) add 
that CMEs can lead to an effectively continuous enhancement 
of planetary radio emission, in particular in young stars with 
their frequent CMEs. 

In spite of the current uncertainties about the rate of exo- 
planetary low-frequency bursts, the possibility of a direct de- 
tection of exoplanets via such bursts is an exciting prospect. 
Being closer to their star than Jupiter is to the Sun, hot Jupiters 
should show much more intense radio en ussion than Jupiter 
itsel f, at least of order 10^ times more ("Lazio et al.', '200 ' 
GrieBmeier et a l.. 2005; Zarka, 2002). In particular, Lazio et al 
( 120041) have computed the expected radio emission for the ex- 
trasolar planets known at the time (mostly hot Jupiters) and 
found that they should have bursts of flux densities close to 
1 mJy at frequencies of several tens of MHz up to a few GHz, 
while a few lower-mass planets should emit at frequencies 
around 0. 1-1 MHz and rea ch tens to hundreds of mJy. Recently, 
GrieBmeier et al. ( 2007bl) have compared predictions for three 
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separate emission processes for currently known exoplanets, 
and find lower peak fr equencies and average emitted powers 



than lLazio et al.1 ( |2004 . 



Referring to Fig. |4] the detection of a 1 mJy bu rst at 5cr sig 



nificance over a typical duration of 15 minutes (ILazio et al 



2004ft at a frequency of 30 MHz and at 25% bandwidth requires 
of order 10^ dipoles, while detecting a 10-mJy burst of equal 
duration at 3 MHz (and again with 25% bandwidth) still re- 
quires about 10'* dipoles; at fixed frequency, the required num- 
ber of dipoles is inversely proportional to the required sensitiv- 
ity. Source variability and their Galactic nature would help to 
overcome the confusion limit. The requirement for the number 
of dipoles is driven by the need to detect a burst within its dura- 
tion. Furthermore, there is no need to map the spatial structure 
of a burst, i.e. the sensitivity requirement is on flux, not surface 
brightness. Hence, the dipoles can be distributed over baselines 
that yield arcminute-scale resolution, which is sufficient for es- 
tablishing positional correspondence to observations at other 
wavelengths. 

Before investing in the construction of such an array, it is im- 
portant to verify that the temporal pulse broadening due to ISM 
turbulence ( ^2. 4. 6b does not smear out an exoplanetary burst 
with a duration of a; 15 minutes beyond recognition. We used 
the Fortran implementation of the NE2001 model for the dis- 
tribution of free electrons in the Galaxy by ICordes and Lazio 
(2002) to determine the scattering regime, pulse-broadening 
times, and angular broadening for 176 exoplanetJ^ The plan- 
ets have a mean distance of 45 pc, and sample a large range 
of Galactic longitudes and latitudes (though their distributions 
naturally follows that of Galactic stars, so that they are concen- 
trated towards the Galactic plane). 

Figure [15] shows the pulse-broadening times and angular 
broadening scales for the exoplanets as a function of distance. 
Although all exoplanets are in the strong scattering regime 
at frequencies below 30 MHz, the angular broadening even at 



' ' Taken from the list at |http : //exoplanets . org| 
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Fig. 15. Pulse broadening times (top panel) and angular broadening scales (bottom panel) cause d by ISM turbul ence towards known extrasolar planets, as 
calculated from the NE2001 model for the distribution of free electrons in the Galaxy (Cordes a nd Lazi 1 l2002h . The temporal and angular broadening is 
shown for the two different frequencies as indicated; the strength of the former scales as y"'* '*, that of the latter as v"--^. The angular broadening is always 
negligible compared to the likely resolution of a lunar array (compare Fig. [5). The temporal broadening at 1 MHz is short compared to the likely exoplanetary 
radio burst duration of k 15 minutes for all but the most distant planets behind those ISM regions with the strongest turbulence. 



1 MHz reaches at most a few 10s of arcseconds, which is neg- 
ligible compared to the likely resolution of a lunar array (com- 
pare Fig. Is). The pulse broadening times at 10 MHz are nearly 
exclusively well below 1 ms, and in all cases less than 10 ms. At 
1 MHz, the pulse broadening is a factor of 10'* '* larger (eq.fTSIl. 
but still below 1 sec and hence negligible for all planets up to 
70 pc distance. Fewer than 1/6 of the known planets are at a 
greater dis tance than 70 pc, a nd only 5 of the 113 planets con- 
sidered by lLazio et a l. are predicted t o emit at frequen- 
cies below IMHz (Griefimei er et al. I l2007bl do not report any 
predictions below 1 MHz because they considered only Earth- 
based radio observatories). Moreover, as noted above, there are 
indications (P. Zarka, priv. comm.) that the pulse broadening 
may increase less rapidly with wavelength than expected from 
eq. ( fT6b . But even if that equation holds, Fig.fTSlshows that an- 
gular and pulse broadening by the ISM impose no significant 
constraints on the observability of exoplanetary radio bursts. 



3.5. Ultra-High Energy particle detection 

3.5.1. Hadronic Cosmic Rays 

Ultra-high energy cosmic rays (UHECRs) are a rare pop- 
ulation of elementary particles, likely composed of pro- 
tons, atomic nuclei, and potentially photons or neutrinos 
dNagano and Watson , I2OOOI) . Their energies can reach much 
higher values than those currently achievable in man-made 



particle accelerators. The total flux of cosmic rays is about 1 
particle nT^ sterad"' year"' above a fiducial particle energy 
of 10'^ eV. The cumulative flux spectrurrP^ drops roughly as 
E^^ with higher energy thresholds, flattening to £ '^^ below 
the so called "knee" in the cosmic ray spectrum at 2 x 10'^ 
eV dHorandel , 2003). Above the knee, the cosmic ray flux is 
roughly given by 

Fcr(E > £c) = 30 particle km-^ sterad"' day"' [-[^f^] (31) 
at£e >2 - lO'^eV. 



At the highest energies, the particle spectrum extends up to 
and perhaps beyond 10^° eV. The currently best information 
on t he spectrum is available from the Pier re Auger Observa- 
tory ( The^iOTe^Aiiger_Conaboration ^ 2008b . showing the well- 
known "ankle" around 4 x lO'^ eV and a steepening beyond 
4 X 10'*^ eV. The latter is likely due to the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff', a suppression of extragalactic UHECRs 
due to proton-gamma interactions with the cosmic microwave 
background. Origin and nature of UHECRs is still largely un- 
known, but a link wi th nearby accreting supermassive blac k 
holes has been made (IThe Pierre Auger Collaborationi l2007h . 
A major obstacle to understanding these most energetic parti- 
cles in the universe is their small flux of less than 1 particle per 



" i.e., the total number of particles above that energy 
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Fig. 16. Radio flux from cosmic-ray or neutrino events interacting w itli the Moon, as observed at 100 MHz with a single dipole located on the Moon. The 
flux and event rate were computed using the formulae given by Scholt en et all j2006h . 



square kilometer per century, requiring large detectors. Here 
the moon, with its large surface exposed to empty space, can 
play a significant role as a particle detector target. 

UHECRs can be detected through their intense particle cas- 
cades which they initiate when encountering a target. For ex- 
ample, if UHECRs hit the lunar surface a particle shower will 
pass through the lunar regolith with a speed greatly exceeding 
the local speed of li ght (the refractive index of the lunar re- 
golith is n = 1 .73, see Alvarez-Muniz et al.l2006|and references 



therein), leading to inte nse Chere nkov light fe as et al 
Daakesamanskii and Zh eleznykhl 1989; Hankin s et al 



1992'; 
1996; 

Gorham et al.L 12004 ). Due to the finite length of the shower 
of a few meters, the emission must be coherent at long radio 
wavelengths. The shower will then show up as a coherent ultra- 
short radio burst on a time scale of At = 1/Av, which is tens of 
nanoseconds up to microseconds in the low-frequency regime 
(v = 1 - 100 MHz). 

Hence, cosmic rays will in principle create an impulsive 
noise background on top of the steady Galactic noise, but 
the emission will in turn also allo w one to them i n mor e de- 
tail. Falcke and Gorhaml(l2003h and lScholten et al. I (l2006h have 
shown that the radio emission from such lunar showers is de- 
tectable at the very highest energies by a LOFAR-like telescope 
even from Earth. Especially for cosmic rays and neutrinos at 
energies > 10^' eV, well beyond the GZK cut-off. Earth-based 
low-frequency observations of lunar radio bursts generated by 
interactions within the moon offer a unique possibility to study 
these elusive particles. 

However, interpretation of these results in the end depends 



on subsurface properties of the moon and shower-development 
at these high energies (e.g., due to the Landau-Pomeranchuk- 
Migdal effect), making in-situ measurements invaluable. This 
is impossible on the Earth since accelerators cannot produce 
particles at such high energies and the Earth's atmosphere pre- 
vents UHECRs from directly hitting the surface. 

Moreover, placing low-frequency antennas on the moon 
would have several interesting advantages compared to Earth- 
bound radio experiments: First, being 400,000 times closer to 
the location of the shower (assuming a fiducial shower dis- 
tance of 1 km on the moon) than the Earth, would, roughly 
speaking, increase the sensitivity to cosmic rays and neutrinos 
by a factor 400, 000^ and allow correspondingly lower-energy 
particles to be detected. Secondly, lunar radio pulses observed 
from the ground get dispersed by free electrons in the Earth 
ionosphere which can be avoided with in-situ measurements. 
Also, the expected radio-quietness at the lunar surface would 
greatly aid in solidly detecting these events. 

Using the equations given bv lScholten et al. I (l2006l) . we find 
that a single dipole antenna operating at 100 MHz with 20 
MHz bandwidth should be able to detect of order 7 cosmic- 
ray events per day at 10'^ eV and above at flux densities of 
100 MJy, corresponding to signal-to-noise levels of 10 or more, 
assuming events at a distance of or less 1 km, arriving from a 
cone of half-opening angle 45° (Fig.[T6b. 

One advantage of placing the radio receivers directly on the 
moon instead of on the Earth is that high-energy events, which 
are rarer, can be detected over greater distances from the loca- 
tion of the antennas. Therefore, a larger fraction of high-energy 
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events is observed, and the differential energy distribution of 
the observed events scales as l/E, even though the intrinsic 
distribution scales as l/E^; the difference arises because more 
energetic events produce quadratically greater fluxes, which are 
detectable over linearly greater distances, increasing the col- 
lecting area quadratically with energy. Thus, a lunar array finds 
a larger fraction of the rarer high-energy events than an Earth- 
based array. 

The total rate of detected events increases logarithmically 
with energy, up to a maximum given by the largest distance out 
to which an event's radio pulse is still visible to the antenna(s). 
The total event count will still be heavily dominated by low- 
energy events, so that the total number of detected events above 
some energy is approximately given by the number of events at 
that energy which occur within the radius for which the signal- 
to-noise ratio S /N meets the detection threshold. 

At least three antennas are necessary to determine the lo- 
cation of a cosmic-ray event by triangulation and hence the 
original energy of the particle. An array of three elements will 
achieve an S/N of about 2.5 times that in Fig. [16] since the S/N 
scales as y/NiN - 1), where is the number of antennas. The 
total event rate above some fixed S /N scales as ln(A^[A^-l]), be- 
cause lower minimum energies E^nin oc 1 / -\/N{N - 1) become 
accessible and the differential event rate scales as 1 

If the observing frequency is lowered from 100 MHz to 10 
MHz, the flux density per event will be much lower, while the 
sky background is slightly higher, and, of course, the available 
bandwidth is smaller All three effects mean that the S/N for a 
fixed event energy is smaller at lower observing frequencies. On 
the other hand, the emission is expected to be more isotropic at 
lower frequencies; while only events above 10^" eV would be 
observable at 10 MHz, all such events are detectable out to a 
radius of 5 km by a single dipole, with a rate of one detectable 
event about every 2 years. This means that observing cosmic 
rays is feasible in the ULW region, and an astronomical array 
of « 100 elements with 5 km spacing could detect of order 30 
events/year at E > 10^** eV. This is in principle competitive 
with the largest cosmic ray detectors today. 

3.5.2. Neutrinos 

The arguments used for deriving the radio flux generated by 
hadronic cosmic -ray events also apply to neutrino interactions, 
so that the basic array parameters are identical for the detection 
of both neutrinos and more massive cosmic rays. The basic dif- 
ference to the hadronic component of UHECRs is that neutrinos 
will penetrate deeper into the moon before they interact. The 

mean free path is 130 km ^Ecr/lO^^eV) and hence the radio 
pulse would come from deep below the surface. Here, lower 
radio frequencies have the additional advantage of penetrating 
deeper into the moon crust, thereb y enlarging the usable detec- 
tor volume. [Scholten et al. (2006) quote an attenuation length 
of 



terial for neutrino interactions. Hence, a simple radio array 
could in principle easily cover several tens of cubic kilometers. 
This would be more difficult at these wavelengths on Earth, 
due the global short-wave radio interference. However, neu- 
trino searches have already been proposed and performed over 
Antarctica at higher frequencies (IBarwick et al l l2006t iHalzenL 
2007; Spiering, 2002). Therefore, the competitiveness of a lu- 
nar array needs to be evaluated more thoroughly. It may also 
be that the difference in geometry of the moon and the density 
of the lunar crust leaves an interesting parameter space open. 

These experiments could run in parallel with an imaging ar- 
ray. They would require an additional signal path with a pulse 
detection algorithm and a raw time series buffer of several mi- 
croseconds to milliseconds per antennas. The trigger event time 
has to be at the tens of nanosecond level. Unlike for most other 
applications, the antennas would require a large acceptance in 
the horizontal planes to detect "grazing-incidence" pulses and 
toward the ground. For simple dipoles without ground plane 
an up-down ambiguity remains; however, bandwidth-limited 
pulses produced in the near field above the lunar surface appear 
exceedingly unlikely, so that this up-down ambiguity should 
not be a concern. The distinction between neutrinos and pro- 
tons/nuclei in the cosmic radiation would be made based on 
the geometry: only neutrinos would come from well below the 
array. 

As a side-effect of a lunar neutrino detector, one could ob- 
serve the Earth as a neutrino source. The emission is due to 
the production of atmospheric neutrinos which are generated in 
particle showers initiated by cosmic rays hitting the Earth's at- 
mospherj^ Since the terrestrial neutrino flux is well known, 
one could wonder about using the Earth-Moon baseline for 
neutrino oscillation measurements. 

3.6. Meteoritic Impacts 

Finally, we point out the possibility that impulsive radio noise 
measurements might also be a potential tool to investigate me- 
teoritic impacts on the moon. Such impacts seem to produce 



optical flashes which are visible even fro m Earth (lOrtiz et al 



2000t ICudnik et all 120031; ICudnikl 1200"/). It is not inconceiv- 
able to think that such impacts are also accompanied by low- 
frequency radio emission. This may be in analogy to strong 
electromagnetic pulses (EMPs), known from nuclear explo- 
sions, which have already been used to explai n radio emissio n 
from cosmic rays in the Earth's atmosphere ( Colgatel 1 1967 ). 
One could also imagine an interaction between the explosion 
and local electrostatic or magnetic field, however, all of this 
remains speculative at the moment. 

At least, radio emission from impacts has already been used 
by radio experiments in free space, where radio antennas ex- 



/ = 9 m (v/GHz)" 



(32) 



If this scaling continues all the way to 10 MHz, then a sin- 
gle dipole can see a volume of ~ 1.5km^ of detector ma- 



" It is amusing to note that the level of neutrino emission from a rocky planet 
like Earth is in fact a signal for the presence of an atmosphere. Showers 
in a solid medium that is not shielded by an atmosphere produce far fewer 
neutrinos due to the much shorter shower length. Unfortunately the level 
of atmospheiic neutrino emission is far too low to be used in the study of 
exoplanets. 
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hibit radio pulses due to high-speed impacts of dust grains. For 
example, based on early findings of Voyager II, a radio antenna 
on the Cassini spacecraft was used to measure dust particles 
from Saturn's rings hitting the entire spacecraft (Kurth et al., 
I2OO6I) . The idea here is that the energy of the impact creates 
an expanding ionized plasma with an ambipolar electric field, 
which in turns creates a strong voltage pulse on the antenna or 
spacecraft. The time scale of the pulse is on the ms scale and 
hence should be well discernible from CR impacts. Moreover, 
the frequency spectrum of the radio pulse seems related to the 
size of the impacting particle. The radio signals might therefore 
reveal the size and mass distribution of the particles. However, 
this effect typically is associated with electrostatic noise gen- 
erated at or near the antenna itself. In how far this effect can 
applied to a lunar antenna remains unclear at present. 



4. Summary & Discussion 

In this paper we have summarized some of the main science 
goals, requirements, and limitations for astronomy at very low 
radio astronomy frequencies below 100 MHz, i.e., ultra-long 
wavelengths (ULW). In this frequency range the moon, in par- 
ticular its far side, offers by far the best (albeit not cheapest) 
site for conducting observations. In several cases a lunar tele- 
scope even offers unique insights. 

The main limitation for ground-based low-frequency tele- 
scopes is degradation of the image quality due to the Earth's 
ionosphere in the 10-100 MHz frequency range. At very low 
frequencies (< 10-30 MHz), the ionosphere is even blocking 
our view of the universe completely. Space and in particular lu- 
nar radio arrays are therefore back in the focus. Here we have 
discussed a number of science areas where this frequency range 
may be of special interest: cosmology with primordial HI at 
z - 10-100, extragalactic surveys of radio galaxies, clusters 
and galaxies, astroparticle physics of cosmic rays and neutrinos 
hitting the moon, transient flares from stars and exoplanets, the 
3D structure of the local interstellar medium, and solar system 
objects. 



4.1. Science goals 

Table [1] summarizes the various science goals and require- 
ments. Clearly of greatest interest is the question of the dark 
ages and HI observations at the highest redshifts. This would 
contain a huge amount of original information about the very 
early universe unaltered by the complicated astrophysical pro- 
cesses associated with cosmic reionization. 

For example, a determination of the exact boundaries of the 
Epoch of Reionization through its global signal could in princi- 
ple be achieved with a single antenna in a dark crater near one 
of the lunar poles. This would provide almost perfect shielding 
from any radio interference and temperature (hence gain) vari- 
ations - ideal conditions for such measurements. However, this 
experiment is also being tried from the ground under somewhat 
less than ideal conditions and the outcome has to be awaited 



(e.g., see the loose upper limit just obtained bv lBowman et al 
l2008l with the EDGES experiment). 

The next step up is the measurement of the HI power spec- 
trum above z > 10 and into the dark ages. Since the HI gas will 
trace directly the dark matter distribution, this is an excellent 
way to measure the primordial density fluctuations which led 
to the formation of large scale structures in the universe. 

Measuring the baryonic acoustic oscillations (BAOs) also 
during the dark ages has several advantages, the most impor- 
tant one being the lack of non-linear processes and biases. 
Non-linear effects, e.g., the rapid growth of ionization bubbles 
from the first stars, affect HI measurements during the epoch 
of reionization, while unknown star formation efficiencies af- 
fect measurements at lower redshifts using galaxies. Hence, 
measurements of BAOs in the dark ages could tackle, for ex- 
ample, the issue of n on-Gaussianity in the density fluctuations 
jCoorav et al. . l2008h . This is a quantity directly revealing in- 
formation about the inflationary phase of the universe. Mea- 
suring the BAOs at z = 50 on a scale of tens of arcmin- 
utes requires of order 10^ antennas, because the amplitude of 
BAOs is at the scale of 10s of microKelvin. If one tries to 
push the resolution further down to arcminute scales in order 
to discriminate between cosm ological models as discussed by 



Loeb and Zaldarriagal (120041) . one again requires at least of or- 



der 10^ antennas at z > 20. This is a somewhat futuristic ex- 
periment for a space array. If this is realized with 5 m crossed 
dipoles one requires about 1000 km of wires, which corre- 
sponds to of order 1 ton of metal for a light weight wire of 
one gram per meter While this is not the limiting factor, the 
power consumption and the electronics might be. Such broad- 
band low-frequency receivers currently require a power of the 
order of several Watts, bringing the total energy consumption 
rate up to several lOOkW plus a correlator - hence this would 
require some serious but not completely outrageous lunar in- 
frastructure to be available. 

The main obstacle for achieving the same science goal from 
the ground is the ionospheric seeing - not because of the resolu- 
tion, but because of calibration problems. The detection of the 
HI signal against the enormous sky background requires a per- 
fect removal of all strong confusing sources and their side lobes. 
Since low-frequency observations are all-sky observations, the 
confusing sources will be seen through all the different patches 
of the ionosphere on the sky, making this a calibration problem 
with many degrees of freedom. Once more, the coming years 
have to show down to which frequencies ground-based arrays 
can calibrate out the obstacles imposed by ionosphere, interfer- 
ence, and foreground. Current experiments do not even attempt 
to push into this high-redshift range, since well-calibrated ob- 
servations at frequencies v < 30 MHz corresponding to 21 -cm 
emission from z > 50 seem completely out of the question for 
ground-based observations. 

As the ultimate goal one could imagine an actual imaging 
experiment of the dark ages, delivering "all one needs to know" 
about the Universe. However, as we have seen here, the required 
collecting area is several tens of square kilometers to achieve the 
necessary surface brightness sensitivity with arcminute-scale 
resolution. This puts such a project in a very distant future. 
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Table 1. Overview of science cases and requirements. 



Topic 


Our§ 






Requirements 










Frequency Resolution & 


Baselines 


Expected signaP 


N(Antennas) 


fexp (5 cr) 






MHz 




km 








Cosmology 
















Global Epoch of Reionization^ 


13.1.11 


50-150 


2n sr 





5-50 mK 


1 


2h-20d 


Global dark-ages signal 


13.1.21 


30-45 


2n sr 





0-40 mK 


1 or more 


0.3-30 yr 


EoR tomography 


13.1.31 


50-150 


^ lO'-l' 


0.7-20 


^ ImK 


10^-10^^ 


yr 


EoR 21 -cm power spectrum 


13.1.41 


50-150 


^ T-T 


0.05-10 


^ 0.3-5 mK 


10^-10^^ 


yr 


Dark-ages 21 -cm power spectrum|3.1.4| 


30-45 


* 20'-2' 


1-20 


^ 0.03-1 mK 


104-10«_^ 


yr 


Extragalactic surveys 


1231 


10 


1' 


0.1-100 


~ 65 mJy^ 


300 


2 yr 


Galactic surveys 


111 














Solar System neighbourhood 


13.3.11 


0.1-10 


degrees 


0.3-30 


lOOOOK 


10-100 


yr 


Origin of cosmic rays 


13.3.21 


0.1-30 


1" 


(3-30)xl03_^ 


155,000K 


100,000 


100 d 


Transients 


m 














Solar and planetary bursts 


13.4.11 


0.1-30 


degrees 


0.5-200 


MJy 


1-100 minutes-hours 


Extrasolar planets 


13.4.21 


0.5-30 


< 1' 


> 35-1000 


1-10 mJy 


104-105_£ 


15 min 


Ultra-High Energy particles 


|33] 


10-100 


N/A^ 


0-5 


100 MJy 


1-100 N/A (bursts)jl 


Meteoritic Impacts 


1221 















" Surface briglitnesses given in Kelvin, fluxes in Jansky 

^ Here we use 10 < z < 20 as fiducial redshift limits for the Epoch of Reionization (EoR), and 30 < z < 50 for the Dark Ages.. 

This requirement is driven by the demanding the necessary surface brightness sensitivity over timescales compatible with the human lifetime, i.e., of order of years. Note the upper limit N(Antennas) 
< ^max = 4.7 X 10^ (d'lVy^ from the shadowing limit (eq. |9}. The airay parameters for 21-cin tomography and power-spectruin observations scale as Wahi y ?cxp S " ' )?^(1 + zy' = const, where S is the desired 
sensitivity. 

Confusion limit from eq. J22> : about 3 X 10* sources are expected to this limit in a half-sky survey. 
Moon diameter: 3,476 km 

' The number of antennas is driven by requiring a S-cr burst detection within the burst typical burst duration of 15 minutes, with ?cxp^Aiit = const. 
^ Near-field triangulation of event sites helps, but is not strictly required. 

*■ Expected event rate for hadronic cosmic rays: ~ 7/day at £ > lO"* eV for 1 dipole at 100 MHz; ~ 30/year at £" > 10^" eV for 100 dipoles with 5 km spacing at lOMHz 



Apart from the fundamental cosmology science goal, other 
topics are also of potential interest. Clearly, Galactic and ex- 
tragalactic surveys are needed to shed light on how the Uni- 
verse looks at these wavelengths. One expects to be particularly 
sensitive to fossil radio sources in clusters, and to be able to 
constrain the total energy output of quasars. As shown by the 
survey equation (l24l i. an extragalactic half-sky survey with ar- 
cminute resolution at 10 MHz is possible within two years with 
an array of a few hundred elements and baselines up to 100 km 
down to a 5 cr detection of the confusion background. 

A fairly unique science application at low frequencies would 
be the 3D tomography of the local warm interstellar medium. 
Since the hot gas in our galaxies starts to become optically 
thick between 0. 1 and 3 MHz, a ULW radio array could see the 
entire spatial distribution of this gas from the Galactic Center to 
the Local Bubble by simply stepping through all frequencies. 
The gas will have imprinted in its structure the recent history of 
the Galaxy and in particular of the solar system environment. 
This is of particular interest since the Local Bubble pressure 
has a direct effect on the size of the heliosphere and thus on 
the cosmic ray flux on Earth and its related effects on weather 
and biological mutation rates (Scherer et al., 2006). Hence, one 
might be able to shed some light on the co-evolution of life in 
the solar system and its interstellar environment. To perform 
this survey, a lunar telescope would require several tens to 
hundred dipoles, several tens of kilometer baselines and a broad 
frequency coverage down to 0. 1 MHz. That could be an realistic 
science goal for the first lunar telescope installations. 

Another new application in ULW astronomy is the search for 
transient sources. Here very little is known a priori. Jupiter is 
the prime example with its cyclotron maser flares at 10 MHz 
which are completely invisible at frequencies above 40 MHz 
due to a sharp spectral cutoff, depressing the flux by 5 orders 
of magnitude (Zarka, 2007). It has been suggested that this 
could be used for the detect i on of the magnetospheres of ex - 
oplanets dLazio et all 12004 iGrieBmeier etall l2005l l2007bh . 
but will certainly also be of interest to investigate solar flares in 
nearby stars. Transient detection requires good instantaneous 
(m, v) coverage, tens of thousands of dipoles, some spatial res- 
olution (baselines of 10-100 km), triggering mechanisms, and 
buffering electronics. 

Due to its huge mass and lack of atmosphere, the moon is 
also a unique ultra-high energy cosmic ray detector. This is 
currently being used by ground-based experiments which ex- 
ploit the fact that particle cascades in the lunar regolith should 
produce strong radio Cherenkov emission. At the very highest 
energies for cosmic ra ys this promises by fa r the most effective 



detection mechanism (^ Scholten et al. . 20061) . The main limit at 



present are systematic uncertainties of lunar properties. These 
could be overcome with the help of a relatively small lunar 
surface radio array which would make in-situ measurements 
of these cosmic-ray-induced radio flashes. Significant advances 
might already be possible with a first prototype array of only 
a few dipoles. Of particular interest is the transparency of the 
lunar regolith below 10 MHz. If this continues to increase with 
decreasing frequency, already a few dipoles could provide an 
effective neutrino detector volume of many cubic kilometers in 



the future. A by-product of such impulsive radio noise measure- 
ments could be the measurement and localization of meteoritic 
impacts on the moon. 



4.2. Open questions 

It is obvious that a full-blown lunar ULW radio array with 
thousands of elements cannot be the first step. Indeed, there 
are a number of open questions about observability and de- 
sign constraints for a moon-based ULW radio telescope, which 
need to be answered before a design for a lunar array can be 
contemplated in detail, let alone be finalized: 

(i) What are the properties of the lunar "ionosphere", in 
particular at the likely sites for a lunar array on the far 
side and near the poles? 

(ii) What are the subsurface and electrical properties of the 
moon? 

(iii) How strong is the shielding of terrestrial radio waves the 
in polar region? 

(iv) How strong it the impulsive radio background noise on 
the surface? 

(v) How strong are electromagnetic effects associated with 
meteorite impacts that could cause impulsive radio noise 
on the lunar surface? 

(vi) Are the IPM and ISM scattering extrapolations from 
higher frequencies still accurate near 3 MHz? 

(vii) What is the detailed synchrotron and free-free emissivity 
and optical depth of the ISM, as a function of galactic 
coordinates? 

(viii) What is the lowest frequency at which ionospheric fluc- 
tuations above a terrestrial observatory can be calibrated 
sufficiently well to achieve the dynamic range necessary 
for imaging faint sources? 

Some of these questions will be answered by LOFAR or MWA 
in the near future, followed later by the SKA. Instruments on 
current-generation lunar orbiters like the on-going Japanese 
KAGUYA/SELENE mission will also provide new insights into 
questions about the electrical properties of the moon. Never- 
theless, some questions will remain that will requke measure- 
ments and prototyping on the moon directly. 



4.3. Stepping stones en route to a large lunar array 

The big advantage of the phased dipole array technology is 
that it is robust, very flexible, and scalable with the number of 
antennas. Therefore, one can easily define a program with sev- 
eral steps signified by a growing number of antennas. Dipole 
antennas can in principle be packaged in a small device and 
therefore could initially piggy-back on other missions for fea- 
sibility studies. The following steps could be envisaged: 

(i) a dedicated pathfinder experiment with a single antenna 
to 

- analyze the transmission of the lunar ionosphere 

- investigate the surface and subsurface electrical prop- 
erties of the Moon 
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- determine the beam shape and stability of a dipole on 
the surface 

- monitor the radio noise environment on the lunar sur- 
face 

- measure the redshift of the global Epoch of Reioniza- 
tion (EEB 

(ii) a two-element interferometer with a moving antenna that 
can 

- demonstrate lunar interferometry 

- perform a first simple ULW all-sky survey 

(iii) a small lunar ULW telescope consisting of > 3 dipoles 
in the 100 MHz regime that can 

- act as cosmic-ray, neutrino, and meteoritic impact de- 
tector (§E5jl[221l 

- demonstrate lunar interferometry 

(iv) an imaging telescope operating in the 10 MHz regime 
with 30-300 dipoles distributed over 30-100 km, for 

- extragalactic surveys of large-scale radio structures in 
clusters, radio galaxies, etc. 03.2b 

- tomographic mapping of the Galactic ISM 03.3. lb 

- monitoring planetary radio bursts ( ^3.4.1b 

- neutrino and cosmic ray measurements ( ^3.5b 

(v) a full-size array optimized for the 10-100 MHz regime, 
with 10^-10^ antennas spread over tens of kilometers to 

- investigate the cosmic dark ages ( ^3. lb 

- search for transient Galactic sources, in particular ex- 
oplanets 03.4.1b 

Among the early prototypes, a simple two-element interfer- 
ometer near one of the poles is a particularly appealing con- 
cept. It could consist of a fixed antenna and one installed on a 
rover, with the rover operating over several months and slowly 
moving further away. The combination of lunar rotation and 
motion of the rover out to tens of kilometers would provide a 
decent (u, v) coverage and allow one to perform an interesting 
imaging survey already with a small mission. 



4.4. Potential realizations 

The realization of these arrays is a matter of current study. 
Dipoles are the preferred primary element for the simplicity 
and ease of deployment, but also magnetic loop antennas and 
tripoles are possible options. Since the Moon's surface is a 
much better insulator than the Earth's, in principle the dipoles 
could be placed dir ectly on the l unar surface, without a defined 
ground plane (Taka hashii l2003h . 

To keep data rates for earth-moon telemetry manageable (de- 
tails below), an array with a large number of elements requires 
a correlator to be placed on t he moon. F o r the d eployment of 
a large number of antennas, IJones et al.l (120071) have consid- 
ered printing of antennas onto very thin foil that can be rolled 
out by astronauts. However, for a smaller number of dipoles 
robotic deployment by rovers or even simple ballistic ejection 
mechanism seem to be sufficient. 

A major concern is the power supply to individual antenna 
elements. This has to be achieved either locally, through small 
solar cells and batteries for some fractional night operation. 



or through wire connections to a central solar power genera- 
tor. The limited capacity of current rechargeable batteries is a 
major limitation for any night operation. Day-time operation 
has always the danger of being affected by strong interference 
from the sun and the lunar ionosphere at the very longest wave- 
lengths. 

A second challenge is the data connection between the in- 
dividual elements and the central data processor ("correlator") 
on the Moon, where the data streams from all antennas have to 
be correlated with each other, and the connection between the 
correlator and the ground station on Earth. For 10 MHz band- 
width and 8 -bit digitization at 20MS/s a connection of x 200 
Mbit/sec is needed. Rather than digitizing the radio waves at 
the dipoles and risking self-generated interference, it may be 
preferable to carry out the digitization close to the correlator. 
For this, the analog radio signal could be modulated onto a 
laser or a glass fibre connection. Due to the curvature of the 
moon and its rocky nature, laser connections will only work 
over a limited distance. Wireless radio communication at high 
frequencies could also be considered, but this option is disfa- 
vored as it would disturb the radio-quiet nature of the site. The 
data rate demands between the correlator and Earth are much 
less than the intra-lunar connections, but have to go through a 
relay satellite for a location on the far side. 

The data correlator itself does not appear to be a funda- 
mental problem. Modern general-purpose F ield Programmable 



Gate Array s (FPGA)-based correlators (e.g., von Herzen','l998] 
Woody et al., 2004; Takeuchi et al., 2005; Parsons et al., 2006 
2007h are already able to integrate the cross-correlation of 30 
antennas with 10 MHz bandwidth essentially on a single chip. 

4.5. Conclusions 

The low frequency end of the electromagnetic spectrum is 
currently not well-explored, and the very-low frequency end is 
indeed one of the last portions of the electromagnetic spectrum 
to remain terra incognita in astrophysics. A lunar low-frequency 
array would therefore enable pioneering observations in a num- 
ber of scientific areas (Tab. [T]i, from the Dark Ages in cosmol- 
ogy to exoplanets and astroparticle physics. The field is still 
in its infancy and serendipity may hold some surprises. After 
all, the transition to long wavelengths also marks a qualitative 
transition, where coherent processes play a much stronger role 
and hitherto unknown ultra-steep spectrum source could reveal 
themselves for the very first time. If lunar exploration becomes 
a reality again in the next decades, then a low-frequency array 
for astrophysics and space science should be an important part 
of it, supplementing the impressive array of new radio tele- 
scopes on the ground. 
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